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ABSTRACT 


Amplification  of  low-li^t-level  optical  images  is  important  for  extending  the  range  of 
lidar  systems  that  image  and  detect  olgects  in  the  atmosphere  and  underwater.  The  use  of 
range-gatmg  to  produce  images  of  particular  range  bias  is  also  important  in  minimizmg 
the  image  degradation  due  to  light  that  is  scattered  backward  fi:om  aerosols,  smoke,  or 
water  along  the  imaging  path.  For  practical  lidar  systems  that  must  be  operated  within 
sight  of  improtected  observers,  eye  safety  is  of  the  utmost  importance.  This  thesis 
describes  a  new  type  of  eye-safe,  range-gated  lidar  sensing  element  based  on  Solid-state 
Raman  Image  Amplification  (SSRIA)  in  a  solid-state  optical  crystal.  SSRIA  can  amplify 
low-level  images  in  the  eye-safe  infi'ared  at  1.556  pm  with  gains  up  to  10^  with  the 
addition  of  only  quantum-limited  noise.  The  high  gains  from  SSRIA  can  compensate  for 
low  quantum  efficiency  detectors  and  can  reduce  the  need  for  detector  cooling.  The 
range-gate  of  SSRIA  is  controlled  by  the  pulsewidth  of  the  pump  laser  and  can  be  as 
short  as  30  -  100  cm,  using  pump  pulses  of  2  -  6.7  nsec  FWHM.  A  rate  equation 
theoretical  model  is  derived  to  help  in  the  design  of  short  pulsed  Raman  lasers.  A 
theoretical  model  for  the  quantum  noise  properties  of  SSRIA  is  presented.  SSRIA  results 
in  higher  SNR  images  throughout  a  broad  range  of  incident  light  levels,  in  contrast  to  the 
increasing  noise  factor  with  reduced  gain  in  image  intensified  CCD’s.  A  theoretical 
framework  for  the  optical  resolution  of  SSRIA  is  presented  and  it  is  shown  that  SSRIA 
can  produce  higher  resolution  than  ICCD’s.  SSRIA  is  also  superior  in  rejecting 
unwanted  sunlight  backgro\md,  further  increasing  image  SNR.  Lastly,  SSRIA  can  be 
combined  with  optical  pre-filtering  to  perform  optical  image  processing  functions  such  as 
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high-pass  filtering  and  automatic  target  detection/recognition.  The  application  of  this 
technology  to  underwater  imaging,  called  Marine  Raman  Image  Amplification  (MARIA) 
is  also  discussed.  MARIA  operates  at  a  wavelength  of  563  nm,  which  passes  most 
efficiently  through  coastal  ocean  waters.  The  imaging  resolution  of  MARIA  in  the 
marine  environment  can  be  superior  to  images  produced  by  laser  line  scan  or  standard 
range-gated  imaging  systems. 


1  Introduction 
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1.1  Motivation  for  Research 

The  desire  to  operate  atmospheric  lidar  (^ht  detection  and  rangmg)  systems  in  the 
“eye-safe”  wavelength  region  around  1.56  pm  has  motivated  much  research  iato  new 
laser  sources  and  the  use  of  nonlinear  optical  elements  to  shift  the  fundamental  laser 
source  to  the  desired  wavelength.  The  region  around  1.56  pm  is  advantageous  because 
this  wavelength  is  not  transmitted  to  the  retina  of  the  eye  as  is  the  1.06  pm  line  of 
Nd:YAG  or  visible  lines  of  other  lasers.  ‘  This  allows  the  use  of  greater  energy  pulses  for 
iUumiaation  of  scenes  of  interest  without  endangering  surrounding  personnel.  Pulsed 
lidar  systems  use  the  technique  of  “range  gating”  to  reject  the  intervening  scattered  light 
due  to  smoke  or  fog.  Range  gating  is  typically  accomplished  by  controlling  the  timing 
between  the  outgoing  laser  pulses  and  the  gating  or  exposure  of  the  sensing  element  of 
the  lidar  system.^ 

The  sensing  element  of  an  imaging  lidar  system  is  usually  an  intensified  CCD  (Charge 
Coupled  Device)  array  or  a  streak  tube  /  CCD  array  combination.^  Image  intensifier 
tubes  operate  by  using  a  photoemissive  detector  coupled  to  a  microchannel  plate  (MCP) 
electron  amplifier.  The  gating  of  the  intensifier  is  accomplished  by  switching  the  voltage 
between  the  photocathode  and  the  input  of  the  MCP  while  the  gain  is  controlled  by 
adjusting  the  high  voltage  across  the  MCP.^  Current  commercial  image  intensifiers  are 
limited  to  gains  of  around  10^,  gating  times  of  around  20  nanoseconds,  and  resolutions  of 
around  50  to  100  pm  at  the  photocathode.  In  addition,  infrared  photocathodes  suffer 
from  low  quantum  efficiencies  that  limit  the  achievable  signal  to  noise  ratio  of  images. 
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The  conteats  of  this  dissertation  introduce  a  new  technique  for  image  amplification 
that  uses  stimulated  Raman  scattering  in  solid-state  crystals.  Solid-state  Raman  Image 
Amplification  (SSRIA)  offers  the  potential  of  extended  gain,  much  shorter  gating  times, 
lower  noise,  and  higher  resolution  than  conventional  image  mtensifiers. 

1.2  Discovery  of  Raman  Scattering 

In  1928,  C.  V.  Raman  discovered  the  inelastic  light  scattering  phenomenon  that  now 
carries  his  name.'*’^  Raman  discovered  that  the  spectra  of  light  scattered  by  liquid 
benzene  included  a  component  that  was  downshifted  in  frequency  from  the  original  light 
source.  He  pointed  out  that  the  frequency  change  was  related  to  the  infrared  vibrational 
frequencies  of  the  molecules  of  the  liquid.  Almost  simultaneously,  the  Russian  physicists 
G.  Landsberg  and  S.  Mandel’shtam  discovered  the  same  effect  in  light  scattering 
experiments  from  aystaUine  quartz.®  For  his  discovery  of  Raman  scattering  and  his  other 
researches  into  the  diffiosion  of  light,  Raman  was  awarded  the  Nobel  Prize  in  1 930 . 

Raman’s  experimental  discovery  was  anticipated  by  theoretical  work  that  predicted 
the  associated  change  of  frequency  in  inelastic  light  scattering.  In  1923,  the  Russian 
physicist,  A.  Smekal,  considered  light  scattering  from  an  energy  conservation  point  of 
view.^  Smekal  showed  that  an  atom  of  mass  Af  in  an  energy  state  Ei  moving  with  a 
velocity  v  can,  upon  collision  with  a  photon  of  energy  Tico ,  make  a  transition  to  another 
energy  state  E],  and  also  experience  a  change  of  velocity  to  v'.  From  the  principle  of 
conservation  of  aiergy  these  quantities  are  related  by  Equation  (1.1) 

—  Mv^  +E,+ha)  =  —  Mv’  +  E,^+fico'  (1.1) 

2  2 


t 
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where  h  is  Planck’s  constant  and  a'  is  the  frequency  of  scattered  light.  Smeckal  showed 
that  the  change  in  velocity  was  negligible  so  that 

E,+fia)  =  E,^+fia>'  (1.2) 


or,  rearranging 


h 


(1.3) 


which  is  interpreted  to  mean  that  the  photon  gives  up  some  of  its  energy  to  the  atom  and 
is  scattered  with  less  energy.  This  is  Raman  scattering.  If  the  change  in  light  frequency 
is  zero  and  the  atom  does  not  change  its  state  I,  then  the  process  is  known  as  Rayleigh 
scattering. 

In  1925,  H.  A.  Kramers  and  W.  Heisenberg  published  a  landmark  paper  in  quantum 
mechanics  that  considerably  amplified  Smekal’s  simple  ideas.®  In  this  paper,  Kramers 
and  Heisenberg  showed  that  not  only  is  it  possible  to  produce  scattered  fight  that  is  red 
shifted  as  described  by  Equation  (1.3),  but  also  the  convCTse  process  is  possible,  whereby 
the  atom  loses  energy  and  the  scattered  fight  is  blue  shifted.  In  modem  terminology, 
Raman  scattering  that  produces  red  shifted  fight  is  called  Stokes  shifted,  for  Sir  George 
Stokes,  who  studied  fluorescence.^  Blue  shifted  fight  is  now  called  anti-Stokes  shifted. 
Kramers  and  Heisenberg  also  introduced  the  concepts  that  fight  scattering  with  frequency 
change  is  incoherent  (in  the  spontaneous  Raman  scattering  regime)  and  that  the  process 
must  occur  by  way  of  a  virtual  intermediate  state  m. 

Raman  made  his  discovery  using  filtered  sunlight  as  the  illumination  source.  After  the 
discovery  and  subsequent  follow-up  papers,  further  work  on  the  Raman  effect  soon 
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transferred  to  the  field  of  chemistry,  where  Raman  scattering  was  exploited  as  an 
analytical  tool.  Progress  in  the  field  seemed  to  stagnate  untU  the  invention  of  the  laser 
brought  a  sudden  revival  of  interest  and  work  in  Raman  scattering.  The  higher  intensities 
and  electric  fields  produced  by  lasers  induce  polarizations  in  a  molecule  or  atom  that  are 
nonlinear  in  the  field.  In  conventional  Raman  scattering  (i.e.  the  type  of  effect  studied  by 
Raman  using  weak  light  sources),  the  strength  of  the  induced  polarization  is  linearly 
proportional  to  the  electric  field  strength.  With  the  higher  field  strength  of  a  laser,  the 
polarization  becomes  nonlinear  in  the  electric  field,  involving  the  square,  cube,  etc.  of  the 
incident  electric  field.  The  most  useful  nonlinear  effect  of  this  type  is  called  stimulated 
Raman  scattering. 

1.3  Stimulated  Raman  Scattering 

Stimulated  Raman  scattering  (SRS)  was  first  observed  in  1962  by  Woodbury  andNg^'^ 
while  studying  Q-switching  of  a  ruby  laser  with  a  nitrobenzene  Kerr  cell.  They  observed 
an  infinred  component  in  the  laser  output  at  a  firequency  that  was  downshifted  by  1345 
cm'^  firom  the  laser  frequency.  Woodbury  and  Eckhardt^^  recognized  that  this  infrared 
component  must  result  from  stimulated  Raman  emission  in  nitrobenzene.  SRS  can  be 
described  classically  as  a  parametric  generation  process  with  one  of  the  electromagnetic 
waves  replaced  by  a  material  excitation  wave  (phonon).  Early  interest  in  SRS  arose 
because  it  could  provide  intense  coherent  radiation  at  new  wavelengths  and  because  it 
was  a  possible  loss  mechanism  in  propagating  high-power  laser  beams  through  the 
atmosphere  or  in  a  laser  fusion  plasma.  Consequently,  most  work  involved  the  study  of 
SRS  in  gases,  both  in  the  atmosphere  and  in  gas  cell  mixtures.  In  particular,  SRS  from 
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H2  gas  has  been  studied  extensively. Typically,  SRS  from  H2  gas  produces  intense 
Raman  output  from  a  gas  cell  of  about  a  meter  long  at  a  few  atmospheres  pressure  using  a . 
visible  laser  intensity  of  several  megawatts  per  square  centimeter. 

1.4  SRS  in  Gas  Cells 

In  the  1980’s,  research  in  the  area  of  laser  fusion  led  to  the  investigation  into  using 
SRS  gas  cells  as  high  power  Raman  amplifiers.  Partanen  and  Shaw‘S  investigated  the  use 
of  low-pressure  molecular  hydrogen  in  a  light  guide  as  a  beam  combiner  for  a  high-power 
KrF  laser.  Such  a  Raman  amplifier  may  be  used  in  the  forward  direction  as  an  optical 
beam  combiner  where  many  pump  beams  are  combined  into  a  single  Stokes  beam.  As 
shown  by  Partanen  and  Shaw  and  others,  such  an  amplifier  need  not  imprint  phase  or 
intensity  variations  onto  the  Stokes  beam,  and  nearly  distortionless  amplification  can  be 
obtained.  Recently,  researchers  at  the  Naval  Research  Laboratory  have  investigated  the 
use  of  a  gas  ceU,  stimulated  Raman  amplifier  for  low-light-level,  quantummoise-limited 
amplification  of  images.*"* 

1.5  SRS  in  Solid-State  Crystals 

Researchers  at  Lite  Cycles,  Inc.,  and  the  University  of  Arizona  Optical  Sciences 

Center  have  been  investigating  SRS  in  solid-state  nonlinear  crystals.  SRS  in  solid-state 

« 

media  offers  the  advantage  of  very  high  Raman  gam  in  a  more  compact,  rugged  design. 
Work  to  date  at  Lite  Cycles  and  the  University  of  Arizona  Optical  Sciences  Center  has 
concentrated  on  the  generation  of  1.56  pm  radiation  through  intracavity  solid-state 
Raman  in  Ba(N03)2  nonlinear  crystals.*^  Intracavity  Raman  lasers  offer  the 


unique  ability  to  convert  a  highly  multi-spatial-mode  pump  beam  into  a  near-dif&action- 
limited  output  beam.  Such  high-brightness,  eye-safe  beams  are  ideal  for  illuminating  the 
object  of  a  pulsed  lidar  imaging  system.  In  addition,  intracavity  solid-state  Raman  lasers 
can  be  designed  to  produce  temporally  short  pulses  through  the  process  known  as 
nonlinear  cavity  dumping.  This  is  also  important  in  lidar  systems,  as  the  range  resolution 
of  a  gated  lidar  system  is  a  function  of  the  convolution  of  the  illuminating  laser  pulse  and 
the  optical  gatewidth  of  the  receiver.  Chapter  2  of  this  dissertation  discusses  the 
dynamics  of  intracavity  Raman  lasers  and  shows  how  short  output  pulses  are  a  natural 
outcome  of  intracavity  Raman  conversion.  A  rate  equation  model  for  the  intracavity 
Raman  laser  is  derived  and  model  calculations  are  presented. 

1.6  Raman  Image  Amplifiers 

The  high  Raman  gain  of  solid-state  nonlinear  crystals  and  the  possibility  of 
distortionless  amplification  suggest  the  use  of  these  crystals  for  the  amplification  of 
optical  images.  Stokes  wavelength  light  reflected  from  an  object  can  be  imaged  through 
a  Raman  crystal  amplifier  that  is  pumped  by  a  pulsed  laser.  The  image  light  acts  as  a 
Stokes  seed  for  SRS  in  the  crystal.  The  initially  weak  image  light  can  be  amplified 
greatly  by  SRS  and  then  relayed  to  an  optical  detector.  In  addition,  because  the  Raman 
gain  in  the  crystal  is  a  function  of  the  pump  intensity,  the  Raman  image  amplifier  can 
produce  range-gated  images  through  control  of  the  timing  between  the  returned  Stokes 
image  puls  e  and  the  amplifier  pump  puls  e. 

Similar  work  has  been  done  using  SRS  gas  cells.  Dimcan^"^  and  others  at  the  Naval 
Research  Laboratory  have  investigated  the  noise  characteristics  and  the  low-light-level 
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imaging  capabilities  of  a  gas  cell,  stimulated  Raman  amplifier.  Their  studies  used 
coUinear  pumping  of  a  l-meter  long  cell  filled  with  molecular  hydrogen  at  20 
atmospheres.  They  used  both  direct  imaging  and  Fourier-transfoim  imaging  of  the 
Stokes  signal  through  the  Raman  amplifier  with  an  amplifia:  gam  of  3  x  lO"^.  Duncan,  et. 
al.,  concluded  that  Raman  amplification  can  be  usied  with  input  signal  levels  approaching 
the  quantum  limit  of  one  photon  per  amplifier  mode.  They  also  concluded  that  Fourier- 
transform  amplification  might  have  advantages  for  minimizing  noise  levels  as  well  as 
potential  advantages  in  maximizing  the  field  of  view. 

Moon,  et.  al.^^,  have  demonstrated  three-dimensional  reflective  image  reconstruction 
through  a  scattering  medium  based  on  time-gated  Raman  amplification  in  a  H2  gas  cell. 
Using  a  firequency  doubled,  Q-switched,  mode-locked  Nd:YAG  laser  producing  30  ps 
pulses  at  532  nm,  they  demonstrated  less  than  1-cm  depth  resolution  through  a  distance 
of  more  than  8  attenuation  lengths  of  turbid  material.  Their  experimental  design  used  an 
optical  delay  arm  to  vary  the  delay  times  of  the  pump  pulse  for  the  Raman  amplifier. 
They  made  the  projection  that  for  an  illumination  level  of  1  J  per  pulse,  lO'*  pixels  per 
image,  and  an  overall  object  reflectivity/collection  efficiency  of  10'^,  an  imaging  depth  of 
~13  attenuation  lengths  should  be  possible.  This  corresponds  to  imagrng  through  27 
meters  of  typical  coastal  water  with  an  attenuation  length  of  2  meters. 

In  Chapter  3  of  this  dissertation,  the  theory  of  Raman  image  amplification  is  presented 
with  a  particular  emphasis  on  solid-state  nonlinear  crystals.  Chapter  3  shows  that  high 
gains  can  be  attained  using  Raman  crystals  of  approximately  1  cm.  x  1  cm.  x  5  cm.  size, 
in  contrast  to  H2  gas  cells  which  require  cells  of  about  1  meter  in  length.  Chapter  3  also 


19 


contains  a  derivation  of  the  noise  properties  of  Raman  image  amplifiers.  Although  the 
Raman  image  amplification  process  does  add  noise  to  the  incoming  signal,  it  will  be 
shown  in  Chapter  3  that  this  noise  is  the  minimum  required  by  quantum  theory  and  is 
more  thau  offset  by  the  high  preamplification  gain  provided  by  SSRIA. 

Chapter  4  contains  a  discussion  of  the  optical  design  aspects  of  Raman  image 
amplifiers  such  as  the  achievable  field-of-view  and  the  limiting  optical  resolution. 

Chapter  5  contains  results  of  an  experimental  demonstration  of  SSRIA.  The 
experiment  consists  of  three  main  optical  parts:  1.)  a  Stokes  seed  source  laser,  which  is  a 
Q-switched,  flashlamp-pumped,  intracavity  Raman  laser  emitting  up  to  30  mJ  per  pulse 
of  1 .56  |rm  radiation  at  a  1  Hz  repetition  rate,  2.)  a  Raman  amplifier  pump  laser,  which  is 
a  Q-switched,  flashlamp-pumped  NdrYAG  laser  emitting  Up  to  89  mJ  per  pulse  of  1.34 
pm  radiation  at  a  1  Hz  repetition  rate,  and  3.)  the  Solid-state  Raman  Image  Amplifier 
which  is  configured  in  a  Fourier  transform  arrangement  where  the  amplifier  crystal  is  at 
the  Fourier  transform  plane  of  the  input  image. 

1.7  Applications  of  Raman  Image  Amplification 

Solid-state  Raman  image  amplification  has  the  potential  to  improve  on  the  current 
state  of  image  intensifier  and  streak  tube  performance,  particularly  in  the  infrared  where 
photocathode  intensifiers  have  poor  quantum  efficiency.  It  offers  the  potential  of  being 
less  expensive,  simpler,  more  rugged,  and  having  higher  performance  than  current  sensor 
devices. 

Applications  for  military  use  include  eye-safe  imaging  lidars  that  can  be  used  for 
autonomous  vehicle  identification  and  targeting  at  20  km  and  greater  ranges,  such  as 
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required  by  the  Air  Force.  Reduced  background  radiation  through  narrow  receiver 
spectral  bandwidths  and  increased  range  resolution  from  short  amplifier  laser  pulsewidths 
are  realizable  with  SSRIA  sensors.  The  increased  range  resolution  will  also  improve 
clutter  rejection.  SSRIA  also  has  the  potential  to  perform  optical  image  processing 
functions  that  will  allow  automatic  target  recognition  or  ahnpoint  designation  in  real-time 
without  the  need  for  extensive  digital  signal  computations.  This  is  readily  accomplished 
in  SSRIA  by  pre-filtering  the  optical  image  before  amplification.  Matched-filtered  output 
from  the  image  amplifier  can  be  obtained.  A  suitable  threshold  can  then  be  applied  to  the 
amplified  image  to  produce  an  all-optical,  automatic  target  recognition  process . 

Other  military  applications  include  underwater  imaging  of  mines  and  submarines. 
Underwater  imaging  of  objects  during  daylight  hours  requires  methods  to  reduce  the 
sunlight  reflected  from  the  water  surface  into  the  camera  or  sensor.  These  reflected 
sunlight  glints  could  overwhelm  a  Hdar  signal  return  from  under  the  water,  even  with 
narrow  bandpass  filters  in  front  of  the  camera  lens.  Sub-nanometer  receiver  spectral 
bandwidths  can  be  realized  with  the  SSRIA  receiver,  enabling  a  system  to  image 
mderwater  objects  with  less  degradation  from  the  sun.  Any  imaging  application  that 
requires  high  range  resolution  -  limited  by  the  amplifier  laser  pulse  width  -  could  be 
realized  from  the  SSRIA  sensor. 

Medical  applications  include  imaging  through  highly  scattering  media.  The  SSRIA 
sensor  would  enable  amplification  and  detection  of  ballistic  photons  transmitted  through 
tissue  with  very  high  range  resolution.  This  resolution  would  enable  production  of  3- 
dimensional  images  of  small,  millimeter-dimension  targets  such  as  tumors.  By  using 
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mode-locking  technology,  the  Raman  illuminator  and  amplifier  pump  laser  pulses  can  be 
temporally  compressed  to  tens  of  picoseconds  FWHM.  The  resulting  gate  for  Raman 
imaging  would  allow  separation  of  the  “ballistic”  photons,  which  pass  unscattered 
through  soft  tissue,  fi-om  the  multiply  scattered  photons,  which  would  normally  mask  the 
desired  image.  Mahon,  et  al.,  have  demonstrated  that  features  as  small  as  125  pm  can  be 
imaged  through  dense  scattering  media  with  light  extinction  factors  of  up  to  e'^^  using 
Raman  image  amplification  in  Hydrogen  gas.  By  using  solid-state  Raman  image 
amplification,  a  medical  imaging  device  can  be  made  smaller  and  more  compact,  with 
even  higher  spatial  resolution.  This  technology  could  result  in  much  earlier  detection  of 
carcinomas,  particularly  those  associated  with  melanomas. 

Chapter  6  discusses  the  optical  design  issues  of  these  applications,  including  issues  in 
the  design  of  eye-safe  imaging  lidar  and  issues  in  the  design  of  an  underwater  lidar 
system  based  on  Raman  image  amplification. 

Chapter  7  presents  conclusions  drawn  from  this  study  and  points  to  future  directions  of 
study. 
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2  Temporal  Dynamics  of  Raman  Lasers 
2.1  Rate  Equation  Model 

The  physical  system  that  is  modeled  in  this  chapter  is  shown  schematically  in  Figure 
2.1-1.  It  consists  of  a  solid-state,  Nd-doped,  YAG  laser  crystal  between  two  mirrors 
which  are  highly  reflecting  at  1.34  pm.  The  Nd:YAG  is  optically  pumped  by  a 
flashlamp.  In  another  optical  cavity  that  shares  one  of  the  highly  reflecting  mirrors  of  the 
1.34  pm  optical  cavity,  there  is  a  Ba(N03)2  nonlinear  Raman  crystal.  The  other  mirror 
of  the  1.56  pm  Raman  cavity  is  a  partially  reflecting  output  coupler.  An  intracavity 
telescope  is  included  in  the  1.34  pm  cavity  that  reduces  the  beam  radius  of  the  1.34  pm 
light  inside  the  Raman  crystal.  A  Pockels  cell  Q-switch  is  also  present  in  the 
fundamental  cavity  to  rapidly  switch  the  optical  loss  in  the  cavity  and  to  produce  giant 
pulse  output. 

HR  @1.34 -1.56  pm 


mirror 
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Consider  first  the  amplification  of  light  in  the  Nd:YAG  rod.  Figure  2.1-2  shows  a 
schematic  diagram  of  the  physical  process  of  light  amplification. 


Incident  radiation 
of  frequency 


Energy  level  diagram 
of  Nd  atom 


Figure  2,1-2  The  propagation  of  radiation  of  frequency  s  Vjj  in  a  medium  of 
atoms  with  a  transition  frequency  . 

The  incident  radiation  is  assumed  to  be  a  quasi-monochromatic  plane  wave  of  frequency 
Vp  and  irradiance  Ip,  {wfcm^)  /  propagating  in  the  +z  direction  (p  represents  “pump” 
since  the  fundamental  at  1.34  pm  is  the  pump  for  the  Raman  process).  The  rate  at  which 
electromagnetic  energy  passes  through  the  plane  of  cross-sectional  area  ^  at  z  is  Ip{z)A . 


At  an  adj acent  plane  at  z  +  Az ,  the  rate  is  Ip{z  +  Az)A  and  the  difference  is : 


[l^{,^Az)-I^{z)]A^j-^(l,A)Az  (2.1) 

in  the  limit  of  very  small  Az .  Equation  (2.1)  gives  the  rate  at  which  electromagnetic 
energy  leaves  the  volume  AAz.  In  terms  of  the  energy  density  Up  of  the  light,  the  rate 


balance  can  be  written: 
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(2.2) 

A  and  Az  are  constant,  and  =/^/c,  where  c  is  the  speed  of  light,  so  Equation  (2.2) 
can  be  written  as: 


This  equation  is  called  the  equation  of  continuity,  and  is  an  example  of  Poynting’s 
theorem,  in  one  space  dimension/^  Equation  (2.3)  applies  to  a  plane  wave  propagating  in 
a  vacuum.  For  a  wave  that  is  propagating  in  a  medium,  the  zero  on  the  right  hand  side  of 
this  equation  must  be  replaced  by  the  rate  per  unit  volume  at  which  electromagnetic 
energy  changes  because  of  the  medium.  This  rate  of  change  is  proportional  to  the  rate  of 
change  of  the  population  inversion  density,  n,  in  the  medium,  and  modifies  the  equation 


of  continuity  to: 


c  dt  dz) 


(2.4) 


where  cr^  is  the  stimulated  emission  cross-section  of  the  pump  {crn^latorn^ ,  and  n  is  the 
population  inversion  density  {atoms j CTr^\ 


In  a  resonant  cavity,  there  are  both  right-  and  left-going  waves  inside  the  medium. 
This  leads  to  two  equations  of  continuity; 


{f^ght  -  going) 

dz  c  at 


(2.5) 
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Adding  these  two  equations  produces: 


If  it  is  assumed  that  there  is  no  vmation  in  ^e  z-direction  (the  so-called  rate  equation 
approximation),  then  the  first  term  in  Equation  (2.7)  can  be  dropped,  resulting  in: 

^(/;+/;)=c<r,,.(';+/;)  (2.8) 

If  the  gain  medium  does  not  fill  the  resonator,  then  n  =  0  outside  the  laser  crystal. 
Integrating  both  sides  of  Equation  (2.8)  over  z  firom  z  =  0  to  z-  (the  length  of  the 
pump  cavity)  we  obtain: 

0  0  ' 

Since  the  LHS  of  diis  equation  is  independent  of  z  (by  the  rate  equation  approxiniation), 
we  obtain  the  relation: 

^(/;+/;)i,=c=r,„;,(/;+/;),  (2.10) 

where  Ip  is  the  length  of  the  gain  medium.  Rearranging  the  terms  in  this  equation  we 
obtain: 
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The  number  of  pump  photons  in  the  cavity,  P,  is  proportional  to  the  total  irradiaace  so, 


nP 


(2.12) 


Relating  the  inversion  density,  n,  to  the  total  nximber  of  inverted  atoms  by  nVg  =  N, 
where  is  the  volume  m  the  gain  medium  which  interacts  with  the  pump  cavity  mode, 
we  obtain; 


dP 

dt 


=  CCJ^ 


(2.13) 


But  Vg  =  lpA,  where  A  is  the  cross-sectional  area  of  the  mode,  and  Vp  =  LpA,  so 
Equation  (2.13)  can  be  written  as; 


dP 

dt 


(2.14) 


Including  the  round-trip  loss  terms  in  the  cavity.  Equation  (2.14)  becomes; 


^  =  l^NP-r^P,  (2.15) 

dt  Vp 

where  Ycp-^pltprt  ^  ^p  is  the  pump  cavity  loss  parameter,  and  tp^,=pjc  is  the 
roundtrip  cavity  time.  The  pump  cavity  loss  parameter  is  given  by; 


<5;,=2a^o^^+ln|^^j  (2.16), 

where  a  is  the  internal  loss  coefficient  for  the  pump,  p^,  is  the  perimeter  of  the  pump 
cavity,  and  R;,  Rz  are  the  reflectivities  of  the  pump  cavity  mirrors. 
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So  far.  Equation  (2.15)  includes  the  effects  of  stimulated  transitions  (stimulated 
emission  and  stimulated  absorption),  and  also  cavity  loss  terms,  but  not  yet  spontaneous 
emission  terms.  The  form  of  Equation  (2.15)  makes  it  simple  to  add  the  effect  of 
spontaneous  emission  to  the  cavity.  A  fundamental  principle  of  quantum  mechanics 
states  that  the  spontaneous  emission  rate  from  any  given  set  of  atoms  into  any  one  cavity 
mode  is  exactly  equal  to  the  stimulated  emission  rate  that  would  be  produced  from  the 
same  atoms  by  one  photon  of  coherent  signal  energy  present  in  the  same  mode.“  So  the 
inclusion  of  spontaneous  emission  results  in  one  “extra  photon”  in  Equation  (2.15); 

(2.17) 

at  Vj, 
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The  corresponding  rate  equation  for  the  atoms  can  be  written  as; 

^  =  (2.18) 

where  is  the  pump  photon  flux  {^photons jcm^  ■  sec) .  The  photon  flux  is  related  to  the 
total  number  of  photons  in  the  cavity  by; 

<D  =—  (2.19) 

p  y 

p 

and  since  nV^^N ,  Equation  (2. 1 8)  becomes; 


The  inclusion  of  the  spontaneous  emission  term  in  Equation  (2.17)  implies  a  loss  of  the 
inversion  density  n.  So  Equation  (2.20)  must  be  modified  to: 
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—  =  -^(P  +  l)iV.  (2.21) 

A  similar  line  of  reasoning  can  be  extended  to  the  Stokes  irradiance  inside  the  Stokes 
cavity.  If  //  and  IJ  are  the  right-  and  left-going  Stokes  irradiances,  then: 

(2.22) 

where  cr^  is  the  stimulated  Raman  scatteriug  cross-section  and  p  is  the  pump  photon 
density  p  =  PjVp  .  Integrating  both  sides  over  z  from  z  =  0  to  (the  length  of  the 

Stokes  cavity),  and  realizing  that  <7^  =  0  outside  of  the  Raman  crystal  of  length  /f ,  we 
obtam: 


(2.23) 


As  done  with  the  pump,  the  number  of  Stokes  photons  in  the  cavity  is  proportional  to  the 
total  intensity,  so  Equation  (2.23)  becomes: 


(2.24) 


Using  the  relation  p  =  PjV^  and  including  the  loss  terms  for  the  Stokes  cavity,  we 
obtain: 


dS  ccT^ 
dt 


PS-rJ 


(2.25) 


P 

For  the  physical  case  considered  in  Figure  2.1-1,  an  intracavity  telescope  is  used  to 
reduce  the  pump  beam  size  in  the  Raman  crystal.  This  requires  a  modification  of 
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Equation  (2.25)  to  account  for  the  increased  pump  beam  irradiance  inside  the  crystal.  As 
was  assumed  previously,  the  total  Stokes  irradiance  is  proportional  to  the  total  number  of 
Stokes  photons  in  the  cavity,  with  the  exact  relationship  given  by; 

s=-^(i:+i:)  (2.26) 

Similarly,  the  expression  for  the  total  number  of  pump  photons  in  the  pump  cavity  is; 

The  effect  of  the  intracavity  telescope  is  to  increase  the  pump  beam  irradiance  in  the 
crystal  by  a  factor  of  Aj  A.^.^  ,  where  4„,.^  is  the  cross-sectional  area  of  the  beam  inside 
the  Raman  crystal.  So  the  pump  photon  number,  P,  in  Equation  (2.25)  must  be  boosted 
by  a  factor  of  AlA.^^^  to  account  for  intracavity  focusing.  The  resulting  expression  is; 


dS  _  ca^  f 

dt  ■  U. 


■  A  "l 

7^  PS-rJ 

,  ‘^inside  ) 


{2.2S) 


Raymer  and  Walmsley“  show  that  the  source  term  for  spontaneous  Raman  scattering 
is  completely  analogous  to  the  case  of  atomic  spontaneous  emission,  i.e.  the  spontaneous 
Raman  emission  rate  is  just  the  stimulated  Raman  scattering  rate  for  one  Stokes  photon  in 


the  cavity  mode.  In  other  words,  we  can  replace  5  in  the  stimulated  emission  term  of 


Equation  (2.28)  with  S+\  to  produce: 


dS  _  c<T^  ^  /,.  Y  ^  ^ 

^dt  Vp  J\Ainside  J 

Neglecting  the  loss  terms  for  the  moment,  the  total  number  of  photons  must  be 
conserved.  This  is  just  a  way  of  enforcing  conservation  of  energy.  Therefore,  the  Stokes 


P(5  +  l)-r.^  (2.29) 
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photon  gain  rate  must  be  balanced  by  a  corresponding  decrease  or  loss  in  the  pump 
photon  rate.  So  the  inclusion  of  pump  conversion  to  Stokes  light  results  in  the  following 
rate  equation  for  the  pump  photons: 

—  =  ^iV(/>+l)-^f-^Y^^V(d'  +  l)-r™P  (2.30) 

dt  ^  j 

Treatment  of  the  second  Stokes  light  in  the  same  fashion  leads  to: 

(2.31) 

where  is  the  stimulated  second  Stokes  Raman  scattering  cross-section.  Integrating 
both  sides  over  z  from  z  =  0  to  z  —  L^  (where  it  is  assumed  that  the  second  Stokes 
radiation  is  contained  in  the  Stokes  cavity),  yields: 

^(/:  +  4) = + 4) .  (2-32) 

which,  when  converted  to  second  Stokes  photons  becomes: 

=  (2.33) 

dt  K  UJ 

Using  the  same  arguments  of  balancing  the  rate  of  loss  of  Stokes  photons,  and  of 
spontaneous  second  Stokes  emission,  leads  to  the  complete  set  of  coupled  rate  equations 
for  the  physical  system  (Equations  (2.21)  and  (2.30)  have  been  repeated  here): 

dt  ^  ^ 


(2.34) 
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dP 

dt 


p  p 


r  1  V 


\^sJ\AnsUkJ 


P[s+\)-r,^p 


(2.35) 


dS  ccr^  1 

f  ^  1 

pf  C  .  l\  ss 

II 

UJ 

V  inside  J 

'  s 

UJ 

'[SS^\)-yJ  (2.36) 


dSS 

dt 


(2.37) 


To  check  the  conservation  of  energy  of  this  system  of  equations,  take  the  sum  of  the 


rates  to  find; 

N  +  P+S +SS  —  —y^pP  —  y  -y 

(2.38) 

or,  equivalently. 

N^-P+S^SS+y,pP  +  yJ  +  y,JS  =  0. 

.  _ 

(2.39) 

This  is  essentially  a  statement  of  the  Manley-Rowe  relations^^  for  this  physical  process. 
If  we  integrate  Equation  (2.39)  over  time,  we  will  find  a  useful  relationship  between  the 
total  energies  of  the  pump,  Stokes,  and  second  Stokes  pulses  and  the  initial  atomic 
inversion  in  the  Nd:YAG  gain  medium; 

N+p+s+ss+rcpP+ VcsS + r  = o  (2.40) 

or, 

[  A(oo)  -  N{0)] + [P(oo)  -  P(0)] + [5(<x>)  -  5(0)] + [55(qo)  -  55(0)]  + 

ao  00  oo 

+  =  p.41) 

0  0  0 


But  the  initial  and  long-time  conditions  are,  respectively: 
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N{0)  =  N.,  P{0)  =  S{0)  =  SS{0)  =  0  (2.42) 

and, 

iV(Qo)  =  P(oo)  =  5(oo)  =  6'^(oo)  =  0  (2.43) 


where  N.  is  the  initial  atomic  inversion.  Since  is  the  number  of  pmnp  photons  per 
second  leaking  out  of  the  system,  {hcjX^y^pP  is  the  power  (in  W)  leaking  from  the 


system,  and  =  —/pj  is  the  total  energy  that  leaked  out  of  the  cavity.  Thus 


Equation  (2.41)  can  be  written  as: 


^  7  7 

(2.44) 

he  he  he 

or, 

2  Ti  hr 

E+^E+^E,,=—N,  (2.45) 

^  X  X  X  ^  ' 

^  p  ^p 

2.2  Numerical  Model 

A  detailed  numerical  model  of  the  intracavity  Raman  laser  dynamics  has  been 
constructed  in  a  Mathematica®  notebook.  Before  discussing  some  of  the  results  of  this 
modeling,  it  is  necessary  to  derive  modified  versions  of  the  coupled  rate  Equations  (2.34) 
-  (2.37).  These  equations  describe  the  rates  of  change  of  numbers  of  inverted  atoms  and 
photons.  The  equations  can  be  converted  to  atomic  and  photon  densities  by  dividing 
through  by  the  appropriate  volumes.  The  relationships  between  atomic  (or  photon) 
numbers  and  atomic  (or  photon)  densities  are  expressed  by: 


N  P  S  SS 

n  =  — ,  p  =  — ,  s  =  — ,  ss  =  — 

K  K  K  K 

S  P  ^  ^ 


(2.46) 


where  n,  p,  s,  and  ss  are  the  atomic,  pump  photon,  Stokes  photon,  and  second  Stokes 
photon  densities,  respectively.  Dividing  through  Equations  (2.34)  -  (2.37)  by  the  proper 
volumes,  as  specified  in  Equations  (2.46)  results  in  the  density  rate  equations; 


where  =  ccr^,  r,  =  ccr^,  and  =  ccr^ .  The  quantity  As  is  the  cross-sectional  area  of 
the  Stokes  mode  volume,  i.e.  A^=VJL^.  The  Stokes  cross-section  is  related  to  the 
running  gain  in  the  Raman  crystal,  gs,  by; 


(2.51) 


The  running  gain  at  a  pump  wavelength  of  1.064  |jm  is  11.0  cm/GW,  and  at  the  pump 


wavelength  of  1.338  pm  it  is  =(1.064/1.338)- ll(cw/GlF)  =  8.75(c/w/GlT).  The 
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Nd:YAG  laser  cross-section  at  1.338  |am  is  <7^  =  1.0  x  lO"'®  .  The  second  Stokes 

cross-section  is  given  by: 


The  set  of  Equations  (2.47)  -  (2.50)  is  solved  numerically  using  Mathematica’s  built-in 
function  for  solving  differential  equations,  NDSolveW?^  The  initial  conditions  are  set  to 
be  the  initial  inversion  density,  and  the  spontaneously  emitted  pump,  Stokes,  and  second 
Stokes  photon  densities.  Table  2-1  shows  a  list  of  the  pump  parameters  used  to  model 
the  intracavity  Stokes  laser  used  as  the  illuminating  laser  (see  Stokes  Laser  Section). 


Table  2-1  Pump  Parameters 


Wavelength 

1.338  |jm 

Cavity  Length 

54.05  cm 

Host  Length 

7.6  cm 

Host  Diameter 

0.6  cm 

Mode  Diameter 

0.6  cm 

Mode  Diameter  in  Raman 

0.24  cm 

HR  Reflectivity 

99.5% 

OC  Reflectivity 

99.5% 

Dichroic  Reflectivity 

95.0% 

Cavity  Lifetime 

5.265  nsec 

Table  2-2  shows  the  list  of  first  Stokes  parameters  used  in  the  model. 


Table  2-2  First  Stokes  Parameters 


Wavelength 

1.556  pm 

Cavity  Length 

20.2  cm 

Raman  Length 

2.2  cm 

Raman  Mode  Diameter 

0.166  cm 

HR  Reflectivity 

99.5% 

OC  Reflectivity 

50.0% 

Cavity  Lifetime 

1.31  nsec 

35 


Table  2-3  shows  the  second  Stokes  parameters  used  in  the  model. 


Table  2-3  Second  Stokes  Parameters 


Wavelength 

1.859  pm 

Cavity  Length 

20.2  cm 

Raman  Length 

2.2  cm 

Raman  Mode  Diameter 

0.166  cm 

HR  Reflectivity 

1.0% 

OC  Reflectivity 

1.0% 

Cavity  Lifetime 

0.14  nsec 

Finally,  shows  the  host  pumping  parameters  used  in  the  model. 


Table  2-4  Host  Pumping  Parameters 


Max.  Flashlamp  Energy 

48.8  J 

Pumping  Efficiency 

1.1  % 

Max.  Extractable  Energy 

59.6  mJ 

Initial  Liversion  Ratio 

2.4 

Final  Inversion  Ratio 

0.3 

Maximum  Intracavity  Irradiance 

2.8  GWW 

Figure  2.2-1  shows  the  model  output  at  the  maximum  flashlamp  energy  of  48.8  Joules. 
As  shown  in  the  figure,  the  total  predicted  Stokes  pulse  output  is  E^{out)  =  33.6  mJ.  The 
Stokes  pulse  shows  an  extremely  fast  rise  time,  and  a  fall  time  which  is  consistent  with 
the  1.31  nanosecond  cavity  lifetime  of  the  first  Stokes  cavity.  Note  that  the  model 
predicts  multiple  pulses.  This  is  due  to  the  fact  that  the  model  calculation  shows  the  first 
Stokes  pulse  reachiag  threshold  and  depleting  the  pump  before  the  pump  pulse  reached 
its  peak  value.  As  a  result,  some  inversion  is  left  over  and  another  pvimp  pulse  builds  up, 
causing  a  second  pulse  of  first  Stokes  light.  This  is  obviously  an  undesirable  phenomena 
for  an  illuminating  laser,  since  the  two  pulses  occur  with  a  time  separation  of  only  about 
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100  nanoseconds.  This  corresponds  to  a  one-way  distance  of  about  15  meters  in  air, 
meaning  that  lidar  return  pulses  from  objects  15  meters  apart  would  coincide  at  the 
receiver,  causing  range-gating  ambiguity. 


0  ,  100  200  300  400 

Time(NanoSec) 


Figure  2.2-1  Model  Output  at  48.8  J  Flashlamp  Input 
Gray  line  -  Undepleted  pump  pulse,  Blue  line  -  Depleted  pump  pulse,  Red  line  - 
First  Stokes  pulse. 

Figure  2.2-2  shows  the  model  output  at  a  flashlamp  energy  of  35  Joules.  This  is  the 
level  at  which  the  seed  Stokes  laser  was  operated  during  the  subsequent  imaging 
experiments.  The  figure  shows  that  only  a  single  first  Stokes  pulse  is  generated  with  an 
output  energy  of  11.1  ml.  The  actual  output  was  measured  to  be  12  mJ  under  these 
conditions  (see  Stokes  Laser  Section).  Thus  the  rate  equation  model  accurately  predicts 
both  the  multiple-pulsing  behavior  and  the  output  energy  of  the  Stokes  pulses. 


Tiine(NanoSec) 

Figure  2.2-2  Model  Output  at  35  J  Flashlamp  Input 
Gray  line  -XJndepleted  pump  pulse,  Blue  line  -  Depleted  pump  pulse,  Red  line 
First  Stokes  pulse. 


3  Raman  Image  Amplifier  Theory 


3.1  Semi-classical  Theory  of  Raman  Amplification 

Raman  scattering  is  a  two-photon  process  whereby  one  photon  at  a)p  is  absorbed  and 
one  photon  at  £0^  is  emitted,  while  the  material  makes  a  transition  from  the  initial  energy 
state,  Eo,  to  the  final  state,  Ei ,  as  illustrated  in  Figure  3.1-1 .  Energy  conservation  requires 
that  -0)^)  =  E^-  Eq=  ticOy ,  which  is  just  the  energy  differaice  between  the  final 
and  initial  states.  A  Raman  active  medium  is  able  to  couple  firequencies  that  differ  by  the 
frequency  of  the  molecular  vibration  or  optical  phonon  in  the  material.  This  coupling 
interaction  between  the  pump  and  Stokes  beams  can  be  analyzed  classically  using  the 
wave  equation  for  the  electric  field  and  treating  the  molecular  vibrations  as  damped 
harmonic  oscillators.*^ 


Figure  3.1-1  Raman  Scattering: 

Pump  photon  absorbed,  Stokes  photon  emitted,  and  material  raised  to 
an  excited  leveL 

By  the  apphcation  of  the  slowly  varying  amplitude  approximation,  the  equations 
coupling  the  pump  and  Stokes  beam  intensities  are; 
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and 


dip 

dz' 


(3.1) 


dz' 


—  S^s^p 


(3.2) 


where  g  is  the  Raman  gain  coefficient.  The  frequency  shifts  and  Raman  gains  for  several 


substances  are  shown  in  Table  3-1 . 


Table  3-1  Frequency  Shift  and  Stimulated  Raman  Gain  for  Several  Substances 


Substance 

Raman  shift  (cm'^) 

Raman  Gain  (cm/GW) 

GasH2 

4155 

1.5(300  K,  10  atm) 

Liquid  N2 

2326.5 

17 

Ba(N03)2 

1047.3 

18.9  (1.56  nm) 

Differential  equations  (3.1)  and  (3.2)  describe  the  depletion  of  the  pump  beam  and  the 


amplification  of  the  Stokes  beam,  respectively,  as  the  two  beams  copropagate.  This  set 
of  equations  can  be  solved  by  considering  that  the  number  of  photons  is  conserved  in 


Raman  scattering.  Thus  the  quantity 

=  (3.3) 


must  remain  constant.  The  solution  to  differential  equations  (3.1)  and  (3.2)  is  given  by 
equation  (3.4); 


4(0) 

4(0) 


exp(gV') 


M$exp(g/oZ') 


(3.4) 


7,(0) 

where  /s(0)  and  /p(0)  are  the  input  inteosities  of  the  Stokes  and  pump  beams, 
respectively.  Figure  3.1-2  shows  the  behavior  of  the  output  Stokes  intensity  as  a  function 
of  input  Stokes  intensity  for  a  fixed  pump  intensity  and  crystal  length.  For  input  Stokes 
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inteasities  that  approach  the  pump  intensity,  the  amplifier  output  Stokes  intensity 
saturates. 


Input  Stokes  Intensity  (W/cm^  ) 


Figure  3.1-2  SSRIA  Output  vs.  Input  Intensity  showing  saturation  behavior 
For  Stokes  beam  intensities  that  are  small  compared  to  the  pump  beam  intensity, 
which  is  the  usual  case  for  the  SSRIA,  equation  (3 .4)  simplifies  to: 

4(^0  =  4(0)  exp(g/pz')  (3.5) 

For  a  SSRIA  crystal  of  length,  Lr,  the  output  Stokes  intensity  is  thus ; 


4(4)  =  4(0)exp(g44)  (^-6) 

The  SSRIA  can  thus  be  characterized  by  a  linear  input-output  relationship, 

(3.7) 

where  the  overall  amplifier  gain,  G,  is  given  by; 


G  =  ExpigIpLj,)  =  Exp{G^). 


(3.8) 
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Figure  3.1-3  shows  tiie  gain  region  investigated  in  this  study.  Typical  amplified 
images  were  made  at  a  gain  level  of  around  900.  The  gain  exponent,  Gr,  is  proportional 
to  the  pump  intensity  in  the  SSRIA  crystal.  By  varying  the  pump  energy,  along  with 
changing  the  focusing  lens  at  the  output  of  the  pump  laser,  the  pump  intensity  in  the 
SSRIA  crystal  was  adjusted  over  the  total  range  shown  in  the  figure.  Note  the  region 
hatched  in  blue  in  the  figmre.  This  represents  a  gain  region  which  is  not  accessible  to 
image  amplification,  because  experiments  showed  that  beyond  a  gain  of  10^  the  amplifier 
begins  to  produce  single-pass  stimulated  Raman  scattering  with  no  signal  input  to  the 
crystal. 


10® 


S  lO'' 

O  900’ 
102 


1 

0  2  4  6  8  10  12  14 

Figure  3.1-3  Available  Gain  Region  of  SSRIA 
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3.2  Minimum  Amplifier  Noise 

One  could  consider  that  an  “ideal”  optical  amplifier  would  be  one  that  amplifies  an 
input  signal  with  the  addition  of  no  extra  noise.  However,  it  has  been  proven  that  it  is 
impossible  to  build  such  an  “ideal”  amplifier.^^  In  fact,  such  a  noise-fi:ee  amplifier  would 
violate  Heisenberg’s  uncertainty  principle,  and  is  therefore  forbidden.  To  show  that  this 
is  trae,  consider  the  usual  statement  of  the  uncertainty  principle  regarding  the 
simultaneous  measurement  of  a  particle’s  momentum,  p,  and  position,  x.  The 
uncertainties,  4p  and  zlx,  associated  with  this  measurement  are  related  by 

ApAx>^  (3.9) 

where  fi  is  Planck’s  constant  divided  by  2n.  If  we  now  consider  photons  propagating 
along  the  x-axis,  the  momentum  of  the  photons  (and  the  associated  uncertainty,  4p) 
given  by 

p=nk=n-^--,tp^—  (3.10) 

c  c  c 

where  kis  the  wave  vector  of  the  electromagnetic  field  associated  with  the  photons,  E  is 
the  photon  energy  and  c  is  tiie  speed  of  light.  The  photon  arrival  time  at  a  particular 
location  (and  the  associated  uncertainty.  At)  are  given  by 

,  =  =  ^  (3.11) 

c  c 

Substituting  the  uncertainties  given  in  Equations  (3.10)  and  (3.11)  into  the  uncertainty 
relation  of  Equation  (3 .9)  yields  the  related  uncertainty  relation  between  energy  and  time: 
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AEAt>-  (3.12) 

2 

If  the  xmcertainty  in  the  energy  is  due  to  an  uncertainty  in  the  number  of  photons  while 
the  photons’  frequency  is  fixed,  then  the  energy  uncertainty  is  AE  =hvAn .  The  phase 
uncertainty  of  the  optical  wave  is  related  to  the  time  uncertainty  by  A^  =  2;rvAt. 
Substitution  of  these  last  two  expressions  into  Equation  (3.12)  yields  the  uncertaiuty 
relation  between  photon  number  and  phase: 

AnA^>—  (3.13) 

Assume  now  that  we  have  an  ideal  linear  optical  amplifier  that  is  completely  noiseless. 
Assume  that  the  photon  nxunber  of  the  input  stream  of  photons  is  ±  A«o .  Then  the 
output  of  the  ideal  amplifier  is 

n±An  =  Gn^±GAnQ  (3-14) 

meaning  that  the  input  signal  is  linearly  amplifier  with  gain  G,  and  no  additional  noise 
has  been  added.  Assuming  that  the  phase  of  the  input  signal  is  jJq  ±  A(Pq  ,  the  output 
phase  is  given  by 

tp±A^  =  ^Q  +  0±Aq>Q  (3.15) 

where  0  is  just  a  phase  shift  accounting  for  the  propagation  through  the  amplifier  and  any 
other  physical  effect.  In  other  words,  this  assumes  that  the  amplifier  is  “phase 
preservmg”.  So  the  net  effect  of  the  amplifier  is  to  increase  the  output  photon  number 
uncertainty  by  a  factor  of  G,  while  leaving  the  output  phase  uncertainty  unchanged. 
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Now  let  us  assinne  that  we  have  an  ideal  optical  detector  that  is  sensitive  to  both 
photon  number  and  phase  and  has  unity  quantum  efficiency.  The  first  condition  implies 
that  the  detector  has  the  smallest  possible  uncertainty,  i.e. 

A«A^  =  l/2  (3.16) 

while  the  second  condition  implies  that  every  photon  that  strikes  the  detector  produces 
one  photocount.^®  Assuming  that  our  ideal  optical  amplifier  is  placed  in  front  of  this 
ideal  detector  to  produce  a  composite,  pre-amplified  photon  detector,  we  can  consider 
this  system  to  be  a  new  input/output  device.  Looking  at  the  input/oulput  relationship  of 
this  composite  device,  we  see  that  the  output  signal  has  the  minimum  possible  uncertainty 
(i.e.  Equation  (3.16)).  To  find  the  corresponding  uncertainty  in  the  input  signal,  we  must 
substitute  Equations  (3.14)  and  (3.15)  into  Equation  (3.16)  to  arrive  at 

This  implies  that  for  arbitrarily  large  gains  G,  the  uncertainty  product  of  the  input  photon 
stream  can  be  made  arbitrarily  small.  This  directly  contradicts  Heisenberg’s  uncertainty 
principle.  Equation  (3.13),  and  so  we  have  reached  a  logical  contradiction  in  our 
argument.  Tracing  back  the  steps  in  the  argument,  we  see  that  the  assumption  of  a 
completely  noise-free  optical  amplifier  must  be  in  error.  Therefore,  such  a  noise-free 
optical  amplifier  physically  cannot  exist  -  some  noise  must  be  added  in  the  amplification 
process  to  satisfy  the  uncertainty  principle. 

This  conclusion  is  valid  for  all  phase-insensitive  linear  optical  amplifiers,  i.e.  linear 
amplifiers  for  which  a  phase  shift  of  the  input  produces  the  same  or  opposite  phase  shift 
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of  the  output.  If  the  phase  shift  of  the  output  is  the  same,  then  the  amplifier  is  called 
phase-preserving-,  if  the  phase  shift  of  the  output  is  opposite,  then  the  amplifier  is  called 
phase-conjugating?^  Examples  of  phase-preserving  amplifiers  are  standard  laser 
amplifiers  and  Raman  amplifiers.  Phase-sensitive  linear  optical  amplifiers  do  not. 
preserve  the  phase  of  the  input  signal  and  have  gains  which  depend  on  the  input  phase. 
An  example  of  a  phase-sensitive  optical  amplifier  is  the  optical  parametric  amplifier 
(OP A).  OPA’s  have  both  a  signal  channel  and  an  idler  channel.  An  OPA  that  has  both 
input  channels  excited  is  sensitive  to  the  phase  of  the  input  signal.  This  phase-sensitive 
amplifier  amplifies  both  the  signal  and  the  noise  in  the  quadrature  that  carries  the  signal, 
but  deamplifies  the  noise  in  the  conjugate  quadrature  that  carries  no  information.^®  In 
this  case,  it  is  possible  to  achieve  “noiseless”  optical  amplification  by  the  concept  of 
“squeezing”.^®  However,  this  can  only  be  accomplished  by  giving  up  the  preservation  of 
the  optical  phase  and  thus  makes  this  method  of  amplification  unsuitable  for  optical 
images. 

Given  that  the  uncertainty  principle  must  always  be  satisfied,  it  is  possible  to  derive 
the  TniTiitnuTn  amount  of  linear  amplifier  noise  consistent  with  it.  This  will  establish  a 
lower  bound  on  amplifier  noise.  Assume  that  we  have  a  linear,  phase-insensitive 
amplifier  which  adds  noise  to  the  signal,  followed  by  an  ideal  detector.  The  tmcertamties 

in  photon  numba-  and  phase  associated  with  the  amplifier  are  An^  and 

respectively.  Likewise,  the  uncertainties  associated  with  the  detector  are  An^  and  . 
The  amplifier  noise  and  the  detector  noise  are  assumed  to  be  uncoirelated,  so  the  total 
uncertainty  of  the  ouqrut  signal  is  given  by 
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An^  =Anl+Anj 
Afp^  =  A(p\^  A(p\ 


(3.18) 


The  uncCTtainty  obtained  at  the  detector  is  related  to  the  input  uncertainty  (before  the 


amplifier)  by 


Ar^  =  G^Ai^ 


Acp^  =  A<pI 


(3.19) 


Substituting  Equation  (3.19)  into  Equation  (3.18)  and  taking  the  product  of  the 
uncertainties  we  obtain 

(3.20) 

By  imposing  the  assumption  that  the  detector  has  minimum  uncertainty,  we  can  substitute 


Equation  (3.16)  into  Equation  (3.20)  to  obtain 


,2.1,  A(p^  f  An^  ^  Anl 


(3.21) 


To  proceed  we  must  make  the  further  requirement  that  the  detector  uncertainty  ratio. 
Anil Aq>\  ,  is  chosen  to  minimize  the  output  number-phase  imcertainty.  This  is  done  by 
minimizing  Equation  (3.21)  with  respect  to  AnljA^l  and  leads  to  the  relation: 


An^  ^ 


(3.22) 


In  other  words,  minimizing  the  output  number-phase  uncertainty  with  respect  to  the 
detector  uncertainty  ratio  means  that  the  amplifier  and  the  detector  need  to  be  matched  in 
terms  of  the  ratio  of  their  introduced  uncertainties.^^  Imposing  the  condition  that  the 


number-phase  uncertainty  product  inferred  for  the  input  signal  is  at  the  minimum 
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allowed,  i.e.  =1/2,  and  substituting  Equation  Q22)  into  Equation  (3.21),  we 

obtain: 

AnlA(pl+An^A(p^-^-^  =  0  (3.23) 

4 

with  the  unique  solution: 

(3.24) 

To  proceed  further,  we  need  to  make  the  assumption  that  the  uncertainty  associated  with 
the  amplifier  is  an  additive  white  noise  process  with  Gaussian  statistics.  If  we  also 
assume  that  the  signal-to-noise  ratio  (SNR)  is  large,  then  the  output  signal  statistics  are 

closely  Gaussian.^^  Based  on  these  assumptions,  the  phase  variance  is  A^l  =  P^llPs , 

where  Ps  and  Pn  are  the  mean  powers  of  the  signal  and  the  noise,  respectively.  The  noise 

>  ' 

power  variance  is  given  by  AP^  =  IPgP^ .  The  noise  power  variance  is  related  to  the 

photon  number  variance  by  AP^  ={hvlff  An] ,  where  T  is  the  detector  sampling 
interval.  Substituting  these  relations  into  Equation  (3.24)  we  obtain 

P«=^j{G-'^)=hvB(G-\).  (3.25) 

where  the  relation  on  the  right  assumes  that  the  sampling  interval  is  J  =  1/25 
(25  =  Bq=  amplifier  FWHM  bandwidth).  This  result  can  be  interpreted  in  the  following 
way:  a5i.»4nTniTig  that  amplifier  noise  is  an  additive  white  noise  process  and  that  the  SNR 
and  gain  G  are  large,  the  miniTTium  possible  amplifier  output  noise  power  corresponds  to 

the  energy  of  half  a  photon  {hv/2)  in  a  given  observation  time  T.  This  interpretation  was 
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also  reached  by  H.  A.  Haus  and  J.  A.  Mullen^®  through  an  approach  which  considered  the 
amplifier  noise  in  the  two  light  field  quadratures  without  any  assumptions  on  their 
statistics.  The  amplifier  output  noise  given  by  Equation  (3.25)  ensures  that  the 
uncertainty  relation  LnLq>  =  1/2  is  satisfied  for  both  the  amplifier  input  signal  and  the 
amplifier  output  signal. 

The  mergy  hvjl  can  be  interpreted  as  being  tbe  zero-point  energy  of  the  quantized 

electromagnetic  field,  which  is  given  by  =  ^+^hv The  minimum  detectable 

noise  power,  given  by  is  often  called  quantum  noise  due  to  its  fundamental 

quantum  origin.  The  (juantum  noise  can  be  contrasted  with  thermal  noise,  which  has  a 

power  of  h^S^|\es.^{hv|kBT)-i^,  where  ks  is  Boltzmann’s  constant  and  T  is  the 
temperature.  At  room  temperature,  thermal  and  quantum  noise  are  equal  at  a  wavelength 
near  44  pm.  At  the  1.556  pm  wavelength  of  interest,  the  quantum  noise  is  13  orders  of 
magnitude  larger  than  the  thermal  noise. 

3.3  Quantum  Theory  of  Noise 

The  quantum  theory  of  noise  in  optical  amplifiers  is  a  sub-field  of  quantum  optics. 
General  information  on  the  field  of  quantum  optics  can  be  found  in  Heitler  ;  Loudon  ; 
Scully  and  Zubairy^^;  Sargent,  ScuUy,  and  Lamb^'*;  and  Meystre  and  Sargent^^.  The 
interaction  of  a  quantized  electromagnetic  field  with  two-level  atoms  can  be  described 
through  the  electric  dipole  interaction  Hamiltonian: 

Kd=eD-E 


(3.26) 
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where  D  and  E  represent  the  atomic  electric  dipole  and  electric  field  operators, 
respectively,  and  e  is  the  electric  charge.  The  quantized  electric  field  is  expressed  in 
terms  of  the  photon  CTeation  and  annihilation  operators  ,  a  through: 

i  =  (3.27) 

where  C  is  a  constant  and  ^  is  a  unit  vector  indicating  the  field  polarization. 

If  the  atomic  quantum  states  are  represented  by  the  kets  |l)  and  l2) ,  corresponding  to 
the  ground  and  excited  states,  with  their  associated  energies,  0  and  hco  ,  then  the  atomic 
electric  dipole  operator  can  be  expressed  in  the  form: 

^  =  i^{|l}(2l+|2}{l|]  (3.28) 

where  dis  a  constant  and  s'  is  a  unit  vector  corresponding  to  the  dipole  orientation.  The 
symbols  |l)(2|  and  |2)(l|  represent  atomic  lowering  and  raising  operators,  respectively. 
Using  Equations  (3.26),  (3.27),  and  (3.28),  the  interaction  Hamiltonian  can  be  written  as: 

+ae-’<“‘‘">llX2|-aV<*''’“>|2Xll}. 

The  interactions  represented  by  each  of  the  four  terms  in  Equation  (3.29)  correspond  to 
absorption  of  a  photon  (terms  proportional  to  a)  or  emission  of  a  photon  (terms 
proportional  to  ),  with  the  associated  changes  in  the  atomic  states. 
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The  electric  field  can  be  represented  by  a  statistical  mixture  of  photon  nmnber  states, 
|«).  The  photon  number  states  correspond  to  the  harmonic  oscillator’s  energy 
eigenstates  and  are  governed  by  the  creation  and  annihilation  relations : 

a^\n)  =  4n  +  \\n  + 

a|«)  =  V«|/i-l) 

A  pure  photon  number  state  has  an  exactly  defined  number  of  photons  but  the  phase  is 
th«i,  necessarily,  indefinite  (see  Equation  (3.13)).  Conventional  laser  sources  do  not 
produce  pure  photon  number  states;  electromagnetic  waves  produced  by  these  sources  do 
not  have  a  definite  number  of  photons.  A  better  quantum  state  representation  for  lasers 
are  the  coherent  states,  defined  as  a  statistical  superposition  of  photon  number  states; 

(3-31) 

With  this  representation,  the  probabihty  of  finding  exactly  n  photons  is  givai  by 

^.=<"|S)  =  kr.  (3.32) 

The  coherent  states  have  a  probability  distribution  which  corresponds  to  Poisson 
statistics. 

In  the  following  simplified  derivation,  we  shall  consider  pure  photon  number  states 
instead  of  a  statistical  superposition.  First,  define  the  atom-photon  number  states 
\i,n){i  =  1,2)  representing  the  whole  atom/field  system.  These  states  verify  Equations 
(3.30)  and  are  orthogonal.  If  we  assume  that  the  initial  state  is  |/,n),  the  probability  of 
finding  the  system  in  the  final  state  \j,m)  after  the  interaction  is : 


(3.30) 
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P.  .  = 

^  jm)in 


(3.33) 

With  aa  initial  state  containing  A:  photons,  the  process  of  photon  absorption  is  then  given 
by  the  probability  process  of  photon  emission  is  given  by  ■ 


Figure  3.3-1  shows  an  energy  diagram  for  the  JBeld  with  four  different  possible 
transitions  associated  with  the  absorption  or  the  emission  of  one  photon  in  a  single 
radiation  mode  of  the  field.  Using  the  definition  for  the  interaction  Hamiltonian, 
Equation  (3.29),  along  with  the  relations  in  (3.30)  and  the  orthogonality  properties  of  the 
atom-photon  number  states,  the  probabilities  associated  with  the  events  shown  in  Figure 


3.3-1  are: 


P 

P 


^nP„ 

=  nP„^i 


(3.34) 


where  the  constants  in  Equations  (3.28)  and  (3.29)  were  ignored,  and  the  atomic  dipole  is 
assumed  to  be  colinear  with  the  field.  The  probabilities  that  are  proportional  to  the  initial 
photon  number  in  the  system,  n,  represent  the  processes  of  absorption  and  stimulated 
emission.  The  +1  term  in  the  last  relation  of  Equation  (3.34)  represents  spontaneous 
emission  because  when  there  are  no  photons  present,  n  =  0;  but  there  is  stUl  a  finite 
probability  of  spontaneously  emitting  a  photon,  i.e.  • 
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FieureSJ-l  Energy  level  diagram  for  the  photon  number  states  |n) 

The  arrows  indicate  the  four  possible  transitions,  wift  the  correspondmg 

prohabiiities,  in  which  one  photon  is  absorbed  or  emitted. 

Now  consider  a  colicction  of  atoms  in  their  groond  and  excited  states  with  density  at 
coordinate  z  of  N,(z)  and  Wz),  respectively,  in  an  infinitesimal  slice  of  medium,  dz.  The 
cross  secUons  for  absorption  and  emission  processes  are  represented  by  <r,  and  tr., 
respectively.  Tire  change  of  the  photon  number  probability  dP.  between  z  and  z*dz  is 

given  by: 

Substituting  in  the  expressions  given  by  Erpadon  (3.34)  and  letting  <t  =  cr,W,  and 
b  =  cr  N,  we  obtain: 

a  I 

^  =  a[nP„.,-{n+l)P„}+b{{n  +  l)P„,,-nP„}  (336) 
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which  is  known  as  the  photon  statistics  master  equation,  and  is  the  starting  point  of  the 
famous  1957  paper  by  K.  Shimoda,  H.  Takahasi,  and  C.  H.  Townes  (STX)  .  This 
equation  is  also  known  as  the  Kolmogorov  equation,  which  governs  a  Markov 

process,  defined  as  the  random  time  evolution  of  a  set  {/^}of  n  variables.  A  Maikov 
process  described  by  forward  or  backward  Kolmogorov  equations  is  also  called  a  birth- 
death-immigration  (BDI)  process. 

The  photon  statistics  master  equation  can  be  used  to  evaluate  the  average  noise  power 
of  an  optical  amplifier.  The  rate  of  change  of  the  average  photon  number  (n)  is  given  by 
the  statistical  mean  of  Equation  (3.36): 

=  =  -n(n+l)P„}+b{Ai(/i+l)P„„ -n"P„}](3.37) 

n  CIZ  CtZ 

Following  STT,  we  assume  that  the  photon  number  is  large  enough  and  use  the 
substitutions  nP„_,^{n  +  l)P„  and  «(n-l)P„,  together  with  the  normalization 

condition  ^  P„  =  1 ,  to  find: 

n 

M~a{{n)+l)-b{n).  (3.38) 

The  solution  of  Equation  (3 .38)  is  then 

=  G(z)(n(0))  +  N{z)  (3.39) 

where  G{z)  is  the  amplifier  gain,  given  by: 

G{z)  =  exp|j[fl!(z')  -b{z')]dz' 


(3.40) 
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and  N[z)  is  the  amplified  noise,  given  by: 

=  p.41) 

0  G{z) 

The  noise  represaited  by  N{z)  comes  firom  the  factor  +1  in  Equation  (3.38)  and  in  the 
photon  statistics  master  equation.  Since  this  factor  was  previously  shown  to  be  attributed 
to  spontaneous  emission,  the  noise  contribution  is  interpreted  as  being  amplified 
spontaneous  emission  or  ASE. 

If  we  assume  that  the  coefficients  a  and  b  are  constant  over  the  amplifier’s  length  (i.e. 
a  uniformly  pumped  amplifier  medium)  and  define 
G  =  G(z)  =  exp[(/2  -  ,  then  Equation  (3.41)  becomes : 

a—b  tjN^  —  N^ 

The  mean  output  noise  power  in  bandwidth  B  can  then  be  given  as: 

'P^=hvxN{z)B=n^hvB{G-^  (3.43) 

where  we  have  used  the  amplifier  spontaneous  emission  factor,  n^ ,  defined  by: 


7^2 

rjN,-N,' 


(3.44) 


At  full  medium  inversion  with  all  of  the  atoms  in  the  excited  state,  N^=0,  and  the 
spontaneous  emission  factor  reaches  its  minimum  value  of  one.  In  this  case  the  mean 
output  noise  power  becomes  =hvB{G~l) ,  which  is  identical  to  Equation  (3.25) 
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found  for  the  miTiiTTmm  noise  allowed  by  quantum  mechanics.  In  other  words,  the 
minimum  amplifier  noise  output  is  obtained  when  complete  population  inversion  is 
achieved  in  the  amplijying  medium. 

In  the  more  general  case  where  the  coefficients  a  and  b  depend  on  the  coordinate  z,  the 
spontaneous  emission  factor  can  be  defined  by: 


=  =  (3.45) 

Gz)-1  Gz)-liGiz') 


G{z)-1  Giz)-liGiz') 

With  this  general  definition,  the  mean  output  noise  power  is  stiU  given  by  Equation(3.43). 


3.4  Photon  Statistics  in  the  Linear  Gain  Regime 

The  photon  statistics  master  equation  has  been  solved  exactly  through  the  use  of  a 
probability  generating  fimction  The  PGFis  defined  as: 

F(x.z)  =  Y.^'-P.{z)  (3.46) 

rt 

and  the  probability  distribution  (z)  can  be  calculated  from  it  by  the  relation; 


nl 


^"F{x,z)' 


dx” 


(3.47) 


The  differential  equation  for  the  amplifier  PGF  is  found  from  Equation  (3.36)  through  the 
expression: 

az  uz 
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where  the  smnmatioiis  m  =  n  +  \  and  m  =  n-\  were  used  to  obtain  a  right  hand  side  that 
is  a  function  of  P„  only.  If  we  use  the  relation  =-xdFldx ,  Equation  (3.48) 
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^F{x,z)'\  ^ 

L 

The  computation  of  these  moments  can  be  found  in  Reference  26,  Appendix  G.  The 
moments  are: 

and 

(if(z))  =  G’(j)[(n’  (0)}  -  (n(0))] + 4G(z)Af(z)(n(o)) + G(z)(/i(0)) (z) + N(z) 

(3.57) 

The  output  photon  variance  is  found  by  using  the  relationcr^(z)  =  ^/7^(z)^-^n(z)^^and 
leads  to  the  important  result: 

cr^(z)  =  G^(z)  cr^(0)-(n(0))  +G(z)^n(0)^  +  iV(z)  +  2G(z)iV(z)(n(o))  +  iV^(z) 

(3.58) 

where  cr^(O)  =  (w^(0))“  (K0))  is  the  input  photon  statistics  variance. 

The  different  terms  that  contribute  to  the  output  variance  of  Equation  (3.58)  are 
typically  called  excess  noise,  shot  noises,  and  heat  noises,  respectively^^.  The  first  term, 
excess  noise,  is  proportional  to  cr^(0)-^n(0)^  and  is  the  only  term  that  depends  on  the 
statistics  of  the  input  light.  If  the  input  signal  follows  Poisson  statistics,  then 
cr^(6)  =  ^n(0)^  and  so  the  excess  noise  is  zero."*”  If,  on  the  other  hand,  the  input  signal 


(«(z)).  (3.55) 


follows  the  statistics  of  chaotic  or  thermal  light,  the  appropriate  probability  distribution  is 
given  by  Bose-Einstein  statistics'*^  with  variance  of  cr^(O)  =  +^A2(0)y  Thus, 
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thennal  light  exhibits  excess  noise  proportional  to  (n(0))^  that  represents  the  classical 
fluctuations  of  the  integrated  intensity.  A  quantity  that  is  commonly  used  to  characterize 
the  amount  of  excess  noise  is  the  Fano  factor,  defined  by: 

/(^)  =  T^-  0-59) 

{n{z)) 

For  Poisson  statistics,  the  Fano  factor  is  /=1,  while  for  Bose-Einstein  statistics, 

/=W+1. 

The  second  noise  contribution,  called  the  shot  noise,  is  equal  to  the  amplifier  mean 
output  power  G(/i(0)^+ iV .  The  term  shot  noise  comes  from  the  effect  of  this  noise  in 
photodetection.  The  shot  effect  causes  photocurrent  noise  with  mean-square  fluctuations 
(i^'^  =  2eB^{i)  where  Be  is  the  electronic  bandwidth.  Since  in  an  ideal  detector, 

(/)  =  ,  where  Bo  is  the  optical  bandwidth,  the  shot  noise  is  proportional  to  the  mean 

number  of  detected  photons.  Actually,  the  use  of  the  term  shot  noise  is  conventional 
here,  since  the  present  analysis  does  not  include  photodetection. 

The  last  two  terms  in  Equation  (3.58),  2G{z)N{z){r^0)'^  +  N^{z) ,  are  called  beat 

noises.  The  first  term,  2G(n(0))A^,  can  be  interpreted  as  the  beat  power  between  the 
electric  field  components  of  the  amplified  signal  and  the  ASE  fallmg  on  a  square-law 
optical  detector  (i.e.  a  detector  which  measures  only  the  intensity  of  the  light).  However, 
this  interpretation  is  unsatisfactory  from  a  quantum  perspective,  as  it  once  again  relies  on 
an  appeal  to  photodetection.  Moreover,  E.  Desurvire^^  shows  that  even  in  the  absence  of 
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spontaneous  emission,  a  noise  term  taking  the  form  appears  at  the  output  of 

an  optical  amplifier.  Thus,  this  term  cannot  be  interpreted  as  a  beat  component,  but  is 
instead  a  result  of  the  statistical  nature  of  the  stimulated  emission  process.  Similarly,  the 
last  term  in  Equation  (3.58),  N^{z),  is  a  result  of  the  statistical  nature  of  the  amplified 
spontaneous  emission  process. 

The  analysis  above  assumed  that  photons  could  only  enter  a  single  mode  of  the 
radiation  field.  For  the  analysis  of  a  Raman  image  amplifier,  this  formalism  must  be 
extended  to  multiple  propagation  modes.  For  this  purpose,  let  M  be  the  number  of 
possible  propagation  modes,  and  let  rt  be  the  total  photon  population  distributed  among 
these  modes.  If  it  is  assumed  that  spontaneous  emission  events  are  equally  likely  for  all 
modes,  the  transition  probabilities  in  Equation  (3.59)  must  be  modified  to: 


Kn.V.2.n={rt  +  n)P„ 


(3.60) 


The  photon  statistics  master  equation  is  then  modified  to^®; 


-^  =  a{(«-l  +  W)P„.,  +  -<}  (3.61) 


and  the  output  mean  and  variance  become: 


(n(r))  =  G(z)(/i(0))  +  W(z)|  (3.62) 

<T^(z)  =  G"  (^)[q-^(0)  -  (w(0))]  +  G(z)(w(0)}  +  m{z) + 2G(z)(«(0))  jV(z)  +  (z) 

(3.63) 

P.  Diament  and  M.  C.  Teich^^  have  examined  the  output  photon  statistics  probability 
distributions  of  optical  amplifiers  with  various  input  signal  statistics.  A  relevant  case  for 
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the  study  of  Raman  image  amplifiers  is  an  input  signal  with  Poisson  statistics.  The  input 
signal  probability  distribution  is  then  given  by"**^: 


-P.(O)  = 


(3.64) 


The  output  distribution  is  found  to  be  described  by  the  noncentral  negative  binomial 
(NNB)  distribution,  also  called  the  Laguerre  distribution; 


N’ 


(N+l) 


n-hfn 


expi- 


^(^o)]  r  (W-1)  [  ^(^o)  1 

iV  +  lJ  V  [  N{N  +  1)\ 


(3.65) 


where  is  the  generalized  Laguerre  polynomial  defined  by'*^: 


{n+a)\ 


j^y!(/i-y)!(y  +  a)! 


{-uy 


(3.66) 


This  result  leads  to  the  conclusion  that  light  amplification  does  not  conserve  the  Poisson 
statistics  for  the  signal.  In  fact,  the  probability  distribution  of  Equation  (3.65)  actually 
describes  the  statistics  of  a  system  where  coherent  and  chaotic  (thermal)  light  fields  are 
superimposed.'*^’'*^  ' 


3.5  Optical  Signal-to-Noise  Ratio  and  Noise  Factor 

SSRIA  is  a  technique  of  gated  optical  preamplification.  There  are  three  main  reasons 
for  preamplifying  optical  signals:  1.)  to  boost  weak  optical  signals  so  that  they  are  above 
the  level  of  thermal  noise  present  in  the  optical  detector  and  its  associated  electronics,  2.) 
to  compensate  for  low  quantum  efficiency  of  the  optical  detector,  and  3.)  to  provide 
optical  gating.  A  requirement  of  preamplification  is  that  it  adds  TninimuTn  noise  to  the 
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incident  signal.  The  noise  performance  of  optical  amplifiers  can  be  characterized  by  the 
power  noweyhctor,  F,  defined  by: 


^  _  SNIC, 


SNR 


(3.67) 


'output 


where  SNR^,is  the  square  of  the  signal-to-noise  ratio  input  to  the  amplifier  and 

SNRj,^.,  is  the  square  of  the  output  signal-to-noise  ratio  from  the  amplifia:.  Here  we  are 
using  the  convention  that  the  optical  signal-to-noise  ratio  is  givai  by: 

SNR  =  ^  (3.68) 

CT 

where  (n)  is  the  mean  optical  signal  in  photon  number  and  cr  is  the  root  mean  square 
(rms)  noise  level,  defined  to  be  the  square  root  of  the  photon  number  variance.  For 
input  signals  with  Poisson  statistics,  the  mean  signal  level  is  (n(0)) ,  while  the  rms  noise 

level  is  (n{Q)'j'^ .  This  leads  to  an  input  SNR  of: 


SNR  = 


m 

V("(0)> 


=^|{4^)) 


(3.69) 


or. 


SNRi^„^={n{0)).  (3.70) 

At  the  output  of  the  Raman  amplifier,  the  square  of  the  SNR  can  be  expressed  as: 


SNRl,^  = 

output 


{{n{L))-{riL))J^  G\L){n{Q)f 


a\L) 


0-^(1) 


(3.71) 
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where  L  is  the  length  of  the  Raman  amplifier  crystal,  represents  the 

mean  photon  numbo-  that  contains  signal  information,  and  =  TllN{z)  is  the 

background  ASE.  Substituting  the  results  of  Equations  (3.62)  and  (3.63)  into  Equation 
(3 .7 1 ),  we  find  that  the  inverse  of  the  output  SNR^  is : 


1  1  fwM  . .  n-2y(i:) ,  w(i)[#(^)+i]' 

W0))f '  ^  0{L)  G"(i)(n(0)) 


(3.72) 


where  /(O)  is  the  input  Fano  factor,  defined  in  Equation  (3.59). 


So  in  the  general  case  of  input  photon  statistics,  the  amplifier  noise  factor  is  found  by 
substituting  Equation  (3.72)  into  Equation  (3.67),  giving: 


p  _ _ 


m 


m-\+ 


\+2N{L)  W(£)[Af(Z.)  +  l]l 

G(L)  *  G=(i)(«(0))  I 


(3.73) 


If  the  input  photon  statistics  are  Poisson  distributed,  the  input  SNR^  is  given  by  Equation 
(3.70)  and  the  Fano  factor  is  unity,  yielding  an  amplifier  noise  factor  of: 


i+2iv  m 

a  *  G\n(0))^  G^{ri0)) 


(3.74) 


If  the  mean  input  is  assumed  to  be  high  enough  that  G(n(0))»  A"  (this  is  equivalent  to 
saying  that  the  mean  input  level  is  much  greater  than  at  high  gain),  then  the  last  two 
terms  drop  out  and  the  noise  factor  becomes: 


^_l+2iV_l+2/i,^((7-l) 

G  G 

In  the  high  gam  limit,  the  optical  amplifier  noise  factor  is  reduced  to: 


(3.75) 
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F(G»1)  =  27>^|  (3.76) 

The  discussion  in  Section  3.3  showed  that  the  limiting  value  of  is  unity,  therefore,  the 

noise  factor  of  a  high  gain  optical  amplifier  is  always  greater  than  2.  If  the  noise  factor 
is  expressed  in  decibels,  then  F>3  dB.  This  is  the  so-called  3  dB  limit  or  quantum  limit 
referred  to  in  the  fiber  communications  literature.'*^ 

To  apply  the  equations  above  to  the  solid-state  Raman  image  amplifier,  we  must  find 
an  expression  for  n^  that  applies  to  Raman  amplification.  Following  the  formalism  of  C. 

Henry  and  R.  Kazarinov'*^,  there  are  two  fundamental  sources  of  quantum  noise  in  a 
system  of  coupled  electrons  and  electromagnetic  fields:  1.)  spontaneous  currents  due  to 
momentum  fluctuations  at  optical  fi’equencies  of  electrons  localized  by  atomic  or 
crystalline  fields,  and  2.)  field  fluctuations  caused  by  the  quantum-mechanical 
uncertainty  of  electric  and  magnetic  fields.  In  quantum  mechanics,  the  current  is  an 
operator  that  is  related  to  the  momentum  operator,  ^ ,  by ; 

J=-P  (3.77) 

m 

where  e  is  the  electron  charge  and  m  is  the  electron  mass  (or  effective  mass).  In  noise 
theory,  the  frequency  spectrum  of  the  fluctuations  of  j  is  given  by  the  Fourier  transform 
of  the  correlation  fijnction  ^y(r)y(0)^ . 

The  other  source  of  quantum  noise  arises  fi'om  the  quantum  uncertainty  of  electric  and 
magnetic  fields.  The  electromagnetic  field  in  a  cavity  can  be  decomposed  into  modes 
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which  can  be  regarded  as  quantum  harmonic  oscillators.  A  mode  has  an  energy  spectrum 
with  a  level  separation  of  h(o  and  a  ground-state  energy  of  hcofl ,  given  by; 

£.=("+jW  (3.78) 

The  fluctuating  electromagnetic  field  in  the  ground  state  referred  to  as  the  field  of 
vacuum  fluctuations. 

Kubo'*^  has  shown  that  the  induced  current  in  an  optical  material  is  related  to  the 
susceptibility  of  the  material  by: 

4®)=^|([y(’').7'(o)]y"'*  (3.79)  . 

where  F  is  a  small  volume  with  dimensions  small  compared  to  the  wavelength  of  light 
but  large  enough  to  include  many  electronic  states.  The  imaginary  part  of  the 
susceptibility  is  given  by: 

=  (3-80) 

— OO 

and  is,  therefore,  the  Fourier  transform  of  the  difference  of  antmormaUy  and  normally 
ordered  correlation  functions  of  spontaneous  current-density  operators  at  two  different 
times.  The  same  Fouriar  integrals  occur  when  we  are  looking  for  the  correlation 

functions  of  the  positive  frequency  components  of  the  spontaneous  current  operator,  , 
which  is  found  by: 


(3.81) 
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This  leads  to  a  statement  of  the  fluctuation-dissipation  theorem  given  by  Callen  and 
Welton***  and  also  derived  by  Landau  and  Lifshitz"^^,  whereby  the  correlation  functions  of 
the  frequency  components  of  the  spontaneous  current  operators  are  separately  related  to 
"  the  imaginary  part  of  the  susceptibility: 

{jJl)  =  ^Z"(g))S((v  -  a)'){n^  + 1) 

(3.82) 

{jl'L)  =  ^X"{o))S[a)-co')ri^ 

where  is  the  average  optical-mode  occupation  factor  of  photons  in  equilibrium  and 
can  be  written  as: 


(3.83) 


where  Nz  and  Nt  are  the  populations  of  the  upper  and  lo  wer  level,  respectively,  and  is 
given  in  Equation  (3.44)  (where  7=1  has  been  assumed). 


The  fluctuation-dissipation  theorem  can  be  generalized  to  correlations  between  the 
operators  at  different  spatial  points."^  In  these  generalized  equations  it  will  be  convenient 
to  express  the  imaginary  part  of  the  susceptibility  in  more  commonly  used  optical 
constants, 

=  =  (3.84) 

where  ^^(x)  is  the  real  propagation  constant,  ^^(x)  is  the  gain  coefficient  (or  the 
attenuation  coefficient),  and  £’"(x)  is  the  imaginary  part  of  the  dielectric  function.  These 


66 


qxiantities  are  related  to  the  real  and  imaginary  parts  of  the  refractive  index  by 
(x)  =  CO n'^{'s)lc  and  (x)  =  2o /z"(x)/c .  The  generalized  fluctuation-dissipation 
theorem  is  then  givei  by: 

j  {Kx)il'(x'))  =  -  l)^(x  - -  ®') 

(3.85) 

(yI'(x')LW)  =  8;r^^„(x)g•„(x)/I,/(x-x')^(^i)-©')• 

The  noise  field  in  a  Raman  amplifier  is  due  to  amplified  spontaneous  Raman 
scattering.  The  optical  frequency  spontaneous  currents  that  are  the  source  of  this  noise 
field  are  due  to  fluctuations  of  atoms  at  optical  phonon  frequencies  mixed  with  the  pump 
field.  This  mixiTig  is  a  result  of  the  nonlinear  interaction  of  the  optical  field  with  die 
atomic  vibrations  in  the  solid.  The  atomic  motion  is  mostly  due  to  the  zero-point 
fluctuations  of  the  vibrational  levels.  Following  the  notation  of  Henry  and  Kazarinov'*^, 
assume  that  each  optical  phonon  mode,  q,  is  executing  simple  harmonic  motion  described 

in  the  usual  way  by  quantized  creation  and  annihilation  operators,  bj  and  b^, 
respectively.  The  energies  of  the  vibrations  are  given  by  the  Hamiltonian, 

p.86) 

where  is  the  angular  frequency  of  mode  q.  The  vibrational  modes  are  in  thermal 
equilibrium,  with  the  operator  averages  given  by: 

(44)=",.  {44)=”, +1.  p-87) 


where 
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(3.88) 


Using  the  interaction  representation,  these  operators  have  a  time  dependence  of: 


4(0  =  4  4^(0  =  (3.89) 

The  Raman  scattering  interaction  can  be  written  as  a  sum  of  products  of  the 

vibrational  mode  amplitudes  the  signal  field  at  the  Stokes  frequency 

+  H.c\ ,  and  the  pump  field  exp(-/o^?) +c.c.] ,  which  is  regarded  as  classical 

and  monochromatic.  Here  H.c.  stands  for  the  Hermitian  conjugate  and  c.c.  stands  for  the 
complex  conjugate.  Calculating  the  spontaneous  and  induced  currents  in  a  small  volume, 
V,  and  keeping  only  resonant  terms  with  products  of  the  fields  that  drive  the  lattice  at  an 
optical  phonon  frequency,  the  interaction  Hamiltonian  can  be  expressed  as : 

=  --h^(t)  +  H.c.  (3.90) 

c 

where  the  coefficients  give  the  strength  of  the  interaction,  and  where  an  effective 

current-density  operator  for  Raman  scattering  has  been  defined.  Using  Equation 
(3.90)  and  Equation  (3 .89),  this  operator  can  be  written  as: 

^  9 


(3.91) 
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The  frequency  components  of  the  Raman  scattering  current-density  operator,  ,  are 
found  by  taking  the  Fourier  transform  of  Equation  (3.91),  see  Equation  (3.81).  The 
correlation  functions  of  are  found  by  using  Equation  (3.89)  and  by  replacing 

with  v\[dNldQ.^)d^^: 


{Jr  CO  Jr  Q)^ 


l7cc^  dN 


V  da 


5{(o-co')n^ 


(3.92) 


Ikc^  dN 
V  da 


The  Raman  gain  is  related  to  the  imaginary  part  of  the  susceptibility,  which  can  be 

written  with  use  of  the  Kubo  formula  (Equation  (3.80))  as  being  proportional  to  the 
difference  in  the  two  correlation  functions  shown  in  Equation  (3.92).  The  result  can  be 
expressed  as: 


gJa,  = 


dN  |2 


(3.93) 


Finally,  the  fluctuation-dissipation  theorem,  which  relates  the  spontaneous  current 
correlation  functions  to  gain,  is  found  by  combining  Equations  (3.92)  and  (3.93)  while 
also  replacing  F"'  by  <5(x  -  x')  : 


f  Atc'^ 

(3.94) 

W)  =  +  l)<^(x  -  x’)s{a)  -  (O') 
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Comparing  tiiis  equation  with  Equation  (3.85),  we  see  that  the  fluctuation-dissipation 
theorem  for  Raman  scattering  agrees  with  the  theorem  describing  a  conventional  laser 
amplifier,  provided  that  we  make  the  substitution  n^p—n^+\.  The  Raman  shift  for 
Ba(N03)2  is  1047.3  cm'\  which  corresponds  to  a  phonon  of  energy 
=2.08x10”^  Joules.  Using  Equation  (3.88),  we  find  that  the  mean  phonon 
occupation  number  at  room  temperature  (300K)  is  =  0.0066 .  Thus,  the  spontaneous 
emission  factor  for  a  Raman  amplifier  is : 

=0.0066  +  1^1  (3.95) 

showing  that  a  Raman  amplifier  has  a  noise  factor  that  is  veiy  close  to  the  ideal  amplifier, 
or,  Fs2.  Thus,  a  solid-state  Raman  image  amplifier  at  high  gain  operates  at  the 
quantum-noise  limit.  Battle,  Swanson,  and  Carlsten^®  came  to  the  same  conclusion  using 
a  different  formalism  incorporating  full  quantum  equations  for  the  Raman  interaction. 

How  does  this  noise  factor  compare  to  the  noise  factor  of  standard  image  intensifiers? 
Assume  the  SNR  of  the  optical  image  input  to  the  Raman  amplifier  is  determined  by  the 
Poisson  noise  of  the  incoming  photons.  Since  for  a  Poisson  random  process  the  variance 
in  photon  number  is  equal  to  the  mean  number  of  photons,  Po,  the  input  SNR  can  be 
written: 

(3.96) 

Substitution  of  the  Raman  amplifier  spontaneous  emission  factor  of  n^=l  into  Equation 
(3.75)  results  in  the  expression  for  the  Raman  amplifier  noise  factor,  Fr: 
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F^=2-— +— =  2-— 
G,  G,  G. 


(3.97) 


where  Gr  is  the  overall  Raman  gain.  If  the  Raman  amplifiCT  is  followed  by  an  optical 
detector  with  quantum  efficiency,  t]  ,  the  overall  amplifier  +  detector  noise  factor  can  be 


written  as: 


51 


Pr  =2 — —+■  ^ 


(3.98) 


Gr  ^Gg 

Thus  the  output  SNR  of  the  combined  Raman  amplifier  and  optical  detector  can  be 
written  as, 


SNR,„„  = 


^output 


1 


Gr  rjG^ 


which  in  the  high  gain  limit  reduces  to: 


(3.99) 


(3.100) 


This  result  shows  that  preamplification  of  optical  signals  by  solid-state  Raman  scattering 
can  compensate  for  low  quantum  efficiency  detectors.  The  net  result  of  preamplification 
followed  by  detection  is  equivalent  to  an  optical  detector  with  a  quantum  efficiency  of 
50%. 


In  the  SWIR  region  region  around  1.56  pm,  photocathodes  are  just  now  becoming 
available  and  their  quantum  efficiencies  are  low  (less  than  20%).  Also,  the  process  of 
electron  multiplication  in  image  intensifiers  results  in  both  image  blurring  and  additional 

r  ■ 

image  noise.  Recent  work  with  intensified  CCDs  has  shown  that  the  noise  factor  of  these 
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devices  is  a  complicated  fimction  of  microcliannel  plate  (MCP)  gain  noise  and  other 
additive  noise/^  The  overall  SNR  performance  of  ICCDs  is  also  degraded  by  the 
quantum  efficiency  of  the  sensing  photocathode.  As  a  result,  the  ouq)ut  SNR  for  an 
ICCD  device  can  be  written  as: 


SNR 


output  ,ICCD 


(3.101) 


rccD 


where  Tjpc  is  the  photocathode  quantum  efficiency  and  Ficcd  is  the  noise  factor  of  the 
ICCD.  Recent  studies  have  shown  that  the  noise  factors  of  current  intensifiers  are  in  the 


range: 


ICCD  = 


1.6 -22for  GEN  II 
35  -  42  for  GEN m 


(3.102) 


and  photocathode  quantum  efficiencies  near  1.56  pm  are  utmost  20%. 

By  taking  the  ratio  of  Equations  (3.99)  and  (3.101),  the  figure-of-merit  (FOM) 
performance  gain  of  the  SSRIA  over  the  ICCD  imaging  system  is: 


FOM  = 


(3.103) 


Giving  the  Raman  amplifier  gain  a  value  of  1000  and  assuming  the  use  of  a  HgCdTe 
array  with  a  quantum  efficiency  of  60%,  equation  (3.103)  gives  an  SNR  performance 
gain  of  at  least: 


FOM  =  2  to  3 


(3.104) 
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over  the  use  of  an  intensifier-CCD  combination.  This  SNR  advantage  directly  translates 
into  a  factor  of  2  to  3  greater  range  for  a  LIDAR  system  viewing  a  resolved  target,  with 
all  other  system  parameters  being  equal. 

Another  advantage  in  using  SSRIA  over  ICCDs  is  the  fact  that  the  noise  factor  of  the 
Raman  image  amplifier  is  not  very  sensitive  to  the  overall  Raman  gain.  In  ICCDs,  the 
noise  factor  of  the  MCP  increases  dramatically  with  lower  gain.  This  limits  the 
TnaxiTniim  SNR  of  an  intensified  image.  Figure  3.5-1  shows  a  comparison  of  noise  factors 
of  the  SSRIA  and  intensified  CCDs  as  a  function  of  gain.  At  lower  gains,  the  noise  factor 
of  SSRIA  is  significant^  lower  than  the  standard  image  tubes. 

Figure  3.5-2  shows  the  overall  SNR  comparison  between  SSRIA  and  ICCDs  as  a 
function  of  input  light  level  (given  in  input  photons  per  detector  pixel).  Tlie  red  fine 
shows  the  upper  limit  of  pure  shot  noise  arising  from  the  Poisson  statistics  of  arriving 
photons.  Shown  in  blue,  the  SNR  of  the  SSRIA  is  just  below  this  ultimate  limit  and 
tracks  the  shot  noise  through  four  decades  of  input  light  level.  The  ICCD  SNR  is  shown 
in  magenta  and  falls  significantly  below  the  SNR  attained  by  SSRIA.  In  addition,  the 
ICCD  SNR  essentially  saturates  at  high  input  light  levels  due  to  the  rising  noise  factor  of 
the  microchannel  plate  (MCP). 


Gain  (G) 


Figure  3.5-1  Noise  Factor  Comparison  of  SSRIA  and  ICCDs 
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Photons  per  Pixel 

Figure  3.5-2  SNR  Comparison  of  SSRIA  and  ICCDs 
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4  Raman  Image  Amplifier  Design 

4.1  Optical  Design  Properties  of  Raman  image  Amplifiers 

Figure  4.1-1  shows  the  optical  layout  for  the  Raman  image  amplifier  used  in  study.  A 
standard  USAF  1951  tti-bar  resolution  target  on  a  glass  slide  was  placed  at  the  focal 
plane  of  the  Fourier  transform  lens.  The  resolution  target  was  trans-illuminated  from  the 
left  by  the  expanded  Stokes  beam.  The  middle  of  the  crystal  was  placed  in  the  back  focal 
plane  of  the  same  lens.  With  this  configuration,  the  two-dimensional  Fourier  transform 
of  the  resolution  target  appears  in  the  middle  of  the  crystal.  The  Former  transform 
distribution  is  then  amplified  by  interaction  with  the  pmnp  beam  in  the  Raman  crystal. 
The  amplified  transform  is  also  in  the- front  focal  plane  of  the  inverse  Fourier  transform 
lens  with  the  IR  vidicon  camera  in  the  back  focal  plane.  Thus  the  inverse  Fourier 
transform  lens  forms  a  real,  amphfied  image  of  the  resolution  target  on  the  camera  plane. 

The  amplifier  pump  beam  was  focused  with  a  +750  mm  lens  and  injected  co-axially 
into  the  crystal  by  a  dichroic  beamsplitter  that  is  a  high  reflector  at  1 .34  pm  and  has  high 
transmission  at  1.56  pm.  At  the  back  end  of  the  crystal,  another  dichroic  beamsplitter 
reflects  the  remaining  1.34  pm  pump  beam  to  a  dump.  In  addition,  several  long-pass 
filters  were  placed  in  front  of  the  camera  to  further  reject  any  remaining  pirtnp  light. 


\  /  Fourier  Transform  ,  Image 

\  /  Plane  Plane 


f=  +150  mm 
Fourier  Transform 
Lens 


Pump  Beam  Pump  Beam 

In  Out 


f=  +150  mm 
Inverse  Fourier 
Transform  Lens 


Figure  4.1-1  Optical  Layout  ofRaman  Image  Amplifier 

Figure  4.1-2  illustrates  the  operation  of  the  SSRIA  in  this  study.  The  resolution  target 


is  called  “image”  in  this  figure.  The  Fourier  transform,  Fr{} ,  of  the  image  appears  in  the 
center  of  the  crystal.  The  Raman  interaction  of  FT  {image}  and  the  injected  pump  beam 
amplifies  FT  {image}  to  Amp[FT{image}],  while  simultaneously  depleting  the  pump 
beam.  The  inverse  Fourier  transform  lens  then  forms  the  amplified  image,  Amp[image], 
at  the  camera  plane. 

There  are  several  things  to  note  about  Figure  4.1-2.  First,  the  pump  beam  distribution 
shown  in  the  figure  is  the  actual  amplifier  pump  beam  distribution  measured  at  the 


Fourier  plane  of  the  SSRIA.  Note  that  the  astigmatism  present  in  the  amplifier  pump 
cavity  results  in  a  non-circularly  symmetric  focused  beam.  The  FT{image}  shows  that 
the  vertical  and  horizontal  resolution  bars  transform  to  vertical  and  horizontal  structure 
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on  the  Fourier  transform  (the  cross-shaped  feature  in  FT {image}).  Higher  spatial 
frequencies,  i.e.  closer  spaced  tri-bars  in  the  resolution  target,  transform  to  correspondmg 
structure  that  is  further  from  the  center  of  the  Fourier  transform.  The  overall  DC  light 
level  of  the  image  is  transformed  to  a  peak  in  the  center  of  the  FT.  As  a  result,  the  pump 
beam  distribution  controls  the  spatial  frequency  content  of  the  jBnal  amplified  image. 
High  spatial  frequencies  that  are  transformed  outside  of  the  pump  beam  distribution 
receive  no  amplification  and  are  lost  in  the  final  amplified  image.  Thus  the  pump  beam 
distribution  in  the  Raman  crystal  is  the  controlling  factor  for  the  Modulation  Transfer 
Function  (MTF)  of  SSRIA.  Note  that  since  the  pump  beam  distribution  is  non-circularly 
symmetric,  the  MTF  in  the  horizontal  and  vertical  directions  will  be  different. 

FT  {Image}  Amp[FT  {Image}] 


Figure  4.1-2  The  Raman  Image  Amplification  Process 
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A  non-circularly  symmetric  pump  beam  is  not  ideal  for  high-fidelity  image 
amplification.  A  solution  to  this  problem  is  illustrated  in  Figure  4.1-3.  By  injecting  the 
pump  at  an  angle  to  a  slab-like  Raman  crystal  and  taking  advantage  of  total  iutemal 
reflection,  the  pump  beam  zig-zags  inside  the  crystal  and  forms  a  more  uniform  gain 
distribution  for  any  input  beam  spatial  distribution. 


Figure  4.1-3  Smooth  Raman  Gain  Distribution  by  Pump  Zig-Zag 
Another  optical  design  consideration  for  S  SRI  As  involves  the  optical  invariant  at  the 

crystal.  The  optical  invariant,  1,  at  a  pupil  plane  in  an  optical  system  is  defined  by: 

I  =  hnu  (1.1) 

where  h  is  the  height  of  the  paraxial  ray,  n  is  the  index  of  refraction,  and  u  is  the  angle  of 
the  chief  ray.  Because  the  center  of  the  Raman  crystal  is  a  Fourier  transform  plane  of  the 
object,  it  is  also  a  pupil  plane  of  the  optical  system.  Figure  4.1-4  shows  a  chief  ray  (in 
red)  passing  through  the  pumped  region  inside  the  Raman  crystal.  The  maximum  angle 
subtended  by  this  chief  ray  is  determined  by  the  geometry  of  the  pumped  region: 


h  2h 


2 


(1.2) 
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since  any  greater  angle  will  result  in  non-overlap  of  the  chief  ray  with  the  pump  beam  (2A 
is  the  diameter  of  the  pumped  region).  Figure  4.1-4  gives  Uie  resulting  expression  for  the 
optical  invariant  at  the  crystal,  where  A  is  the  cross-sectional  area  of  the  pumped  region. 
Since  the  optical  invariant  must  be  the  same  everywhere  in  the  optical  system,  the  optical 
invariant  at  the  entrance  pupil  of  a  receiving  telescope  coupled  to  the  SSRIA  must  be 
equal  to  the  invariant  at  the  crystal.  The  expression  for  the  optical  invariant  at  the 
receiving  telescope  is  shown  in  Figure  4.1-5.  It  shows  that  the  invariant  at  the  telescope 
is  the  product  of  the  half  field-of-view  angle  and  the  radius  of  the  telescope  optic.  Thus, 
the  pumping  geometry  in  the  SSRIA  crystal  is  directly  connected  to  the  FOV  and 
aperture  size  that  can  be  used  in  a  lidar  system. 


Figure  4.1-4  Optical  Invariant  at  the  SSRIA  Crystal 
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For  example,  if  we  would  like  to  image  tanks  and  other  vehicles  at  a  range  of  2 
kilometers,  the  total  field  imaged  by  the  system  should  be  20  meters  to  encompass  the 
entire  vehicle.  Thehalf  FOV  is  then  107n/2000/w  =  0.005ra<f .  A  reasonable  aperture  size 
for  this  demonstrator  is  a  telescope  diameter  of  30  cm.  So  the  optical  invariant  is 
I  =  0.015cm  and  the  required  energy  per  pulse  for  an  image  amplification  of  1000  is  then 
Ep  =  143/71/. 

4.2  Optical  Resolution  of  Raman  Amplifiers 

Figure  4.2-1  gives  a  summary  of  findings  on  the  optical  resolution  of  the  SSRIA  and 
compares  SSRIA  with  the  intensified  CCD.  In  most  ICCD  systems,  the  resolution  is 
limited  by  the  combmed  effects  of  input  and/or  output  fiber  optic  bundles,  the  electron 
optics  of  the  input  and  output  MCP  spacings,  and  microchannel  plate  (MCP)  blur. 
Current  Gen  n  wafer  tubes  have  limiting  resolutions  of  25-32  line  pairs  per  millimeter. 

The  SSRIA  optical  resolution  is  governed  by  the  size  of  the  pumped  region  in  the 
Raman  amplifier  crystal.  This  is  due  to  the  Fourier  transform  relationship  between  the 
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optical  image  plane  and  the  craiter  of  the  crystal  plane.  Optical  resolution  is  t5^ically 
specified  by  the  Modulation  Transfer  Function  (MTF)  which  is  defined  by: 


MTF  = 


•^max  •^min 


(1.3) 


where  Imca  and  Imin  are  the  maximum  and  muiimum  irradiance  values,  respectively,  of  the 
image  of  a  sine-wave  pattern.  The  limiting  resolution  is  usually  defined  as  the  spatial 
firequency  at  which  the  MTF  is  3%.  This  is  considered  the  miniTnum  contrast  which  can 
be  seen  by  the  eye.  The  advantage  in  using  MTF  to  characterize  resolution  is  that  the 
MTF  for  each  separate  part  of  an  optical  system  can  be  found  separately  and  then  the 
product  of  these  separate  MTFs  gives  the  MTF  of  the  entire  system. 


SSRIA 

Resolution  is  Camera-Limited 
Eye-safe  systems  require 
camera  sensitive  at  1 .5  p^m 

-  V idicon  and  IR-CCD  ’  s  have 
resolution  of  ~  ^  Jtm/pixel 

Fourier-plane  imaging  pump 
beam  profile  is  the  effective 
MTF 

-  Crystal  dimension  is 
effective  MTF  when  wave- 
guided  pump  geometry  is 
utilized 


ICCD 

Resolution  is  MCP-Limited 

Electron  optics  in  MCP 
defines  the  ultimate  resolution 

-  Gen  n  ICCD  blur 
size  is  ~  125  p m/pixel 


Figure  4.2-1  Resolution  Comparison  of  SSRIA  and  ICCD 
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The  MTF  contribution  of  the  Raman  Image  Amplifier  can  be  found  by  a  consideration 
of  the  Fourier  transform  relationship  between  the  object  plane  and  the  center  of  the 
crystal  plane.  Figure  4.2-2  shows  two  1 -dimensional  sine  wave  objects  on  the  left  and 
their  corresponding  Fourier  transforms  on  the  right. 


6  cycles  per  mm 

◄ - ^ 


Fourier 

Transform 


3  tiith  pumped 
region 
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100 

8o! 

50 

»»oj 

0! 


-1.5  0  L5mm 


3  mm  pumped 
region 


12.8  cycles  per  mm 


Figure  4.2-2  Fourier  Transforms  of  Sine  Wave  Objects 
The  first  object  has  a  spatial  frequency  of  6  cycles  per  mm  and  the  second  object  has  a 
higher  spatial  frequency  of  12.8  cycles  per  mm.  The  Fourier  transforms  are  seen  to 
consist  of  three  distinct  peaks:  a  central  peak  that  corresponds  to  the  DC  light  level,  a 


positive  frequency  peak  on  the  right,  and  a  negative  frequency  peak  on  the  left.  The  two 
outer  peaks  thus  contain  the  high  frequency  information.  As  seen  from  the  bottom 
transform  pairs,  higher  spatial  frequencies  are  mapped  further  from  the  center  of  the 
Raman  crystal  until,  at  the  cutoff  frequency  of  the  SSRIA,  the  outer  peaks  are  at  the  edge 
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of  the  pump  beam  in  the  mystal.  Higher  spatial  frequencies  than  the  cutoff  frequency  are 
mapped  outside  of  the  pump  region  and  are  not  amplified.  The  cutoff  frequency  of  the 
SSRIA  is  controlled  by  the  diameter  of  the  pumped  region,  d,  and  the  focal  length  of  the 
inverse  Fourier  transform  lens,/2,  by  the  equation: 


P  SSRIA  ~ 


(1.4) 


4.3  Sunlight  Background  Rejection 

Figure  4.3-1  shows  a  schematic  comparison  of  the  spectral  characteristics  of  SSRIA 
and  ICCDs.  For  a  typical  ICCD  imaging  lidar  system,  the  use  of  a  photocathode  implies 
a  broad  spectral  response.  To  limit  the  backgrormd  light  due  to  sunlight,  a  narrow 
bandpass  optical  filter  is  used  in  front  of  the  ICCD  optics.  However,  because  of  the  shift 
of  the  center  passing  wavelength  with  angle  inherent  in  multi-layer  dielectric  bani^ass 
filters,  in  order  to  work  with  a  reasonable  field-of-view  the  bandpass  must  be  limited  to 
around  10  nanometers  FWHM.  Any  narrower  and  the  FOV  must  be  decreased 
accordingly.  Figure  4.3-1  shows  that  the  background  light  within  a  10  nanometer 
bandpass  is  amplifiedfry  the  ICCD  and  appears  as  noise  in  the  image. 

The  spectral  width  of  the  Raman  amplification  process  is  on  the  order  of  0.1 
nanometer  and  is  not  sensitive  to  angular  shift.  Thus,  the  sunlight  background  is  reduced 
by  a  factor  of  1 00  when  using  SSRIA  compared  to  ICCDs. 
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5  Raman  Image  Amplifier  Results 

The  most  significant  findings  and  results  produced  by  this  research  are; 

SSRIA  is  an  eye-safe  LIDAR  technology  operating  at  1 .56  pm. 

SSRIA  can  amplify  low-level  IR  images  by  a  factor  up  to  10^. 

SSRIA  adds  only  quantum-limited  noise  to  the  amplified  image. 

SSRIA  can  compensate  for  low  quantum  efficiency  IR  detectors  and  reduces  the  need  to 
cool  detectors  to  reduce  thermal  noise. 

SSRIA  provides  range-gated  LIDAR  images  with  temporal  rangergates  as  short  as  2  -  10 
nanoseconds. 

A  LIDAR  system  using  SSRIA  will  result  in  an  image  signal-to-noise  ratio  (SNR)  gain  of 
more  than  2  to  3  over  current  intensified  CCD  technology. 

Optical  resolution  of  SSRIA  can  be  much  higher  tiian  ICCDs. 

SSRIA  reduces  the  amount  of  sunlight  background  in  gated  images  by  a  factor  of  10  or 
more. 

SSRIA  can  be  used  to  perform  optical  image  processing  functions  such  as  image 
enhancement  and  automatic  target  detection/recognition. 
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5.1  Stokes  Illuminator  Laser 

5.1.1  Stokes  Illuminator  Laser  Design 

The  Stokes  seed  laser  resonator  design  was  based  on  high  peak-,  high  average-power, 
eye-safe  Raman  resonators  previously  demonstrated.  Using  data  obtained  from  existing 
high  power,  eye-safe  systems,  and  applying  well-known  scaling  laws,  the  multimode 
spot-size  required  in  the  Raman  crystal  was  sized  to  obtaia  the  best  energy  conversion 
efficiency. 

Considering  these  results,  the  optimized  resonator  was  designed  using  the  generalized 
resonator  design  code,  called  Lasica®.  The  resonator  design  optimization  included  (1) 
maximizing  the  energy  extraction  efficiency,  (2)  minimiziTig  the  accrued  round-trip 
aberrations,  (3)  minimizing  the  potential  for  optical  damage,  (4)  mode-matching  the 
Raman  cavity  to  the  pump  cavity,  and  (5)  selection  of  the  proper  Raman  output  coupling 
percentage  -  Raman  gain  combination  to  ensure  the  desired  temporal  pulse  profile. 

Figure  5.1. 1-1  shows  the  layout  of  the  Stokes  laser.  It  is  a  Q-switched,  intracavity 
Raman  laser  that  uses  a  nonlinear  optical  crystal  of  Ba(N03)2  to  shift  the  output 
wavelength  from  the  fundamental  1.34  pm  line  of  NdrYAG  to  the  eye-safe,  Stokes- 
shifted  wavelength  of  1.56  pm.  The  Nd:YAG  rod  is  pumped  by  a  flashlamp  and  the 
fundamental  stable  resonator  cavity  at  1.34  pm  is  defined  by  a  1  meter  radius-of- 
curvature  (ROC)  high  reflectivity  mirror  at  the  rear  end  and  a  10  meter  ROC  high 
reflectivity  mirror  at  the  front  end.  This  10  meter  ROC  mirror  is  shared  by  the  Stokes 
resonator  cavity.  The  1.34  pm  light  path  is  shown  in  red  in  Figure  5.1.1-1.  The  Stokes 
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resonator  cavity,  which  has  a  1.56  |jm  light  path  shown  in  blue,  consists  of  the  10  meter 
HR,  the  crystal,  and  the  flat,  50%  reflective,  output  coupler.  A  dichroic  beamsplitter  at  a 
30°  incidence  angle  couples  the  two  cavities. 


HR  @1.34 -1.56  pm 
10  meter  ROC 


All  Dimensions  in  mm 
(unless  otherwise  stated) 

Figure  5.1.1-1  Stokes  Laser  Layout 

Figure  5. 1.1 -2  shows  the  single-mode  beam  radius  in  the  fundamental  cavity  as  a 
function  of  distance  from  the  rear  HR.  Note  that  the  intracavity  telescope  acts  to  deorease 
the  beam  radius  inside  the  nonlinear  crystal  and  thus  mcreases  the  intensity  of  the  pump 
beam  in  the  crystal.  The  limiting  aperture  in  the  fundamental  cavity  is  the  rod  itself  and 
so  at  the  higher  pumping  levels  (up  to  48.8  Joules  lamp  energy),  the  fundamental  beam  is 
highly  multi-mode  and  extracts  maximum  energy  from  the  rod.  Figure  5. 1.1 -3  shows  the 
multi-mode  beam  radius  as  a  function  of  distance  from  the  reai'  HR  (dashed  red  line). 
Also  shown  in  the  figure  is  the  single-mode  beam  radius  for  comparison  (solid  red  line). 
Note  that  the  beam  radius  in  the  Raman  crystal  is  1 .2  mm. 
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Intracavity 

Ba(N03)2 

Q-S^tch  Nd:YAG  Rod 

Telescope 

f - ^ 

Crystal 

Figure  5.1.1-2  Single-mode  Beam  Size  of  134  ^im  Resonator 


Intracavity  Ba(N03)2 

Q-Switch  Nd:YAG  Rod  Telescope  Crystal 


Figure  5.1.1-3  Multi-mode  Beam  Size  of  1.34  jlm  Resonator 
An  important  consideration  in  the  design  of  these  intracavity  Raman  lasers  is  that  the 
multi-mode  spot  size  inside  the  crystal  must  provide  a  sufficiently  high  intensity  to  reach 
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the  Ramaa  laser  threshold  and  should  be  matched  in  size  to  the  siugle-mode  Stokes  beam 
size.  Due  to  the  phenomenon  of  Raman  beam  clean-up,  the  Stokes  beam  produced  by 
intracavity  Raman  lasers  is  typically  close  to  the  fundamental  Gaussian  mode  of  the 
Stokes  cavity  regardless  of  the  mode  structure  of  the  pump  beam.  Thus  for  maximum 
energy  coupling  and  efficiency,  the  multi-mode  pump  radius  should  match  the  single¬ 
mode  Stokes  radius.  However,  this  is  not  always  possible  because  of  the  coupled  cavity 
dynamics.  Too  much  pmnp  intensity  inside  the  crystal  can  cause  the  creation  of  multiple 
Raman  pulses  or  of  Raman  cascading  to  the  second  Stokes  wavelength.  Figure  5. 1.1 -4 
shows  the  Stokes  single-mode  beam  radius  inside  the  Stokes  cavity.  The  single-mode 
beam  radius  is  0.83  mm,  leading  to  an  area  mismatch  efficiency  of 

=f^T=0.48  (5.1) 

Ypump  y  V  / 

With  the  use  of  an  unstable  resonator  for  the  Stokes  cavity,  this  efficiency  can  be  made  to 
approach  one.  For  this  study  the  efficiency  was  judged  adequate. 


89 


Ba(N03)2 

Crystal 


Cavity  Axis  (nun) 

Figure  5.1.1-4  Stokes  Cavity  Design  (Single-mode) 

5.1.2  Stokes  Laser  Characterization 

The  Stokes  seed  laser  was  built  to  verify  and  characterize  the  design.  The  output 
power  and  temporal  pulse  profiles  were  measured  and  monitored  for  stability.  Figure 
5. 1.2-1  shows  the  spatial  beam  profiles  of  both  the  fundamental  beam  at  1.34  pm  and  the 
Stokes  beam  at  1 .56  pm.  The  pump  beam  shows  the  expected  highly  multi-mode  spatial 
structure.  The  Stokes  beam  is  very  close  to  the  single  transverse  mode  or  TEMoo  spatial 
structure.  The  pump/Stokes  area  mismatch  is  clearly  seen  in  this  figure  and  matches  the 
size  predicted  by  the  spatial  model. 
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Pump  Stokes 


Figure  5.1.2-1  Pump  and  Stokes  Beam  Profiles 
Figure  5. 1.2-2  shows  the  measured  input/output  characteristic  curve  for  the  seed 
Stokes  laser.  The  output  energy  per  pulse  was  measured  at  a  pulse  repetition  frequency 
of  1  Hz.  At  the  maximum  lamp  input  energy  of  48.8  Joules,  the  output  energy  at  1.56  pm 
was  30  milliJoules.  For  the  SSRIA  experiments  described  in  Section  2.5,  the  lamp  input 
energy  was  set  at  35  J,  corresponding  to  an  output  Stokes  energy  of  12  mJ. 


Input  Lamp  Energy  (J) 

Figure  5.1.2-2  Stokes  Laser  Input/Output  Characteristic 

Figure  5. 1.2-3  shows  the  temporal  pulse  profiles  of  the  fundamental  and  Stokes 
beams.  The  fundamental  pump  beam  rises  rmtil  the  Stokes  conversion  threshold  is 
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achieved.  At  that  point  the  pump  is  rapidly  converted  to  the  Stokes  wavelength  and  the 
Stokes  pulse  rises  very  sharply.  The  Stokes  pulse  fuU  width  at  half  maximum  (FWHWQ 
is  seen  to  be  6.9  nsec.  The  Stokes  pulse  depletes  most,  but  not  all,  of  the  population 
inversion  in  the  NdiYAG.  This  is  evinced  by  the  second  hump  of  the  pump  beam. 


5.2  Raman  Amplifier  Pump  Laser 

5.2.1  Raman  Amplifier  Pump  Laser  Design 
The  Raman  amplifier  pump  laser  design  is  a  flashlamp-pumped  Nd:YAG  power 
oscillator  configuration  to  produce  the  pump  power  at  1.34  |mi.  The  resonator  design 
was  accomplished  using  the  generalized  resonator  design  code,  Lasica®.  The 
oscillator/amplifier  design  included  (1)  maximizing  the  energy  extraction  efficiency,  (2) 
minimizing  the  accrued  oscillator  round-trip  and  amplifier  single-pass  aberrations,  (3) 
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rrdnimizmg  the  potential  for  optical  damage,  (4)  matching  the  cross-sectional  area  of  the 
SSRIA  crystal  to  ensure  optimum  coupling  and  (5)  suppression  of  parasitic  1.06  pm 
oscillation. 

Figure  5.2. 1-1  shows  the  layout  of  the  amplifier  pump  laser.  It  is  a  Q-switched 
Nd: YAG  laser  operating  at  the  fimdamental  wavelength  of  1 .34  pm.  The  Nd: YAG  rod  is 
pumped  by  a  flashlamp  and  the  fundamental  stable  resonator  cavity  at  1.34  pm  is  defined 
by  a  20  meta-  radius-of-cm-vature  (ROC)  high  reflectivity  mirror  at  the  rear  end  and  a  4 
meter  ROC,  50%  reflectivity  output  coupler  at  the  front  end.  Also  included  in  the  cavity 
are  a  Pockels  cell  Q-switch,  4  YAG  plate  polarizers,  and  two  mirrors  that  are  highly 
reflective  at  1.06  pm  and  highly  transmissive  at  1.34  pm.  These  two  mirrors  are 
necessary  to  suppress  laser  oscillation  at  1 .06  pm. 


HR  @  1.34  -  1.56  pm 
20  meter  ROC 

Pockels 
CeU 

Q-Switch 


4  YAG  Plates 
3  mm  Thick 


HR  @  1.06  pm 
HT  @  1.34  pm 


Output  Coupler 
R  =  50%  @1.34  pm 
4  meter  ROC 


Nd:YAG  Rod 
6  mm  Diam. 
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^  r  -xi. 

's  I  .  -N^ 
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45  45  50 
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All  Dimensions  in  mm 
(unless  otherwise  stated) 

Figure  5.2.1-1  Amplifier  Pump  Laser  Layout 
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Figure  5. 2. 1-2  shows  the  smgle-mode  beam  radius  in  the  pump  laser  cavity.  The  red 
curve  shows  the  beam  radius  in  the  x-axis  (i.e.  in  the  horizontal  plane,  parallel  to  the 
optical  table)  and  the  blue  curve  shows  the  beam  radius  in  the  orthogonal  axis.  The 
introduction  of  the  YAG  plate  polarizers  and  1 .06  |rm  mirrors  cause  a  small  amoimt  of 
astigmatism  in  the  optical  cavity. 


Q-Switch  Nd:YAG  Rod 


Cavity  Axis  (mm) 

Figure  5.2.1-2  Pump  Laser  Single-mode  Cavity 
Figure  5. 2. 1-3  shows  the  multi-mode  cavity  plot  for  the  pump  laser.  The  limiting 
aperture  in  the  cavity  is  the  rod  itself  and  so  at  the  higho’  pumping  levels  (up  to  70  Joules 
lamp  energy),  the  fundamental  beam  is  highly  multi-mode  and  extracts  maximum  energy 
from  the  rod. 
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Q-Switch  Nd:YAGRod 


Cavity  Aids  (mm) 

Figure  5.2.1-3  Pump  Laser  Multi-mode  Cavity 
5.2.2  Pump  Laser  Characterization 

The  designed  Raman  amplifier  pump  laser  was  built  to  verify  and  characterize  the 
design.  The  output  power  and  temporal  pulse  profiles  were  measured  and  monitored  for 
stability.  Figure  5.2.2-1  shows  the  spatial  beam  profile  of  the  amplifier  pump  beam  at 
1 .34  pm.  The  ptunp  beam  shows  the  expected  highly  multi-mode  spatial  structure. 


Figure  5.2.2-1  Pump  Laser  Beam  Profile 
Figure  5. 2.2-2  shows  the  measured  input/output  characteristic  curve  for  the  amplifier 
pump  laser.  The  output  energy  per  pulse  was  measured  at  a  pulse  repetition  frequency  of 
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1  Hz.  At  the  maximuin  lamp  input  aiergy  of  70.0  Joules,  the  output  energy  at  1.34  pm 
was  95  milliJoules.  However,  as  is  evident  from  the  change  of  slope  in  the  figure,  at 
around  45  J  lamp  energy  some  parasitic  oscillation  at  1.06  pm  starts  to  rob  energy  from 
the  1.34  pm  beam.  For  this  reason,  most  of  the  experiments  described  in  Section  5.3 
used  43  J  of  input  lamp  energy  resulting  in  less  than  80  ml  of  pump  energy  per  pulse. 


Figure  52.1-1  Pump  Laser  Input/Output  Characteristic 


Figure  5. 2.2-3  shows  the  temporal  pulse  profile  of  the  amplifier  pump  beam  as  taken 
by  the  oscilloscope  (the  vertical  scale  is  inverted  -  increasing  pulse  intensity  is 
downward).  The  pump  beam  is  shown  in  green  and  the  Stokes  pulse  is  shown  in  red  on 
the  same  time  scale.  The  measured  FWHM  of  the  amplifier  pump  pulse  is  shown  to  be 
43.15  nsec  on  the  right  side  of  the  figure.  Note  that  the  horizontal  scale  is  10  nsec  per 


division. 
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(Grem  -  amplifier  pump  pulse;  Red  -  Stokes  pulse) 


5.3  Raman  image  Amplification  Resuits 

5.3.1  SSRIA  Design 

Figure  5.3. 1-1  shows  the  SSRIA  Demonstrator  layout  that  was  used  to  demonstrate  the 
feasibility  of  SSRIA.  The  Demonstrator  consists  of  three  main  optical  parts:  1)  the 
Stokes  seed  source  laser,  which  is  a  Q-switched,  flashlamp-pumped,  intracavity  Raman 
laser  emitting  up  to  30  mJ  per  pulse  of  1.56  pm  radiation  at  a  1  Hz  repetition  rate,  2)  the 
Raman  ampUfier  pump  laser,  which  is  a  Q-switched,  flashlamp-pumped  NdrYAG  laser 
emitting  up  to  89  mJ  per  pulse  of  1.34  |jm  radiation  at  a  1  Hz  repetition  rate,  and  3)  the 
Solid-state  Raman  Image  Amplifier  which  is  configured  m  a  Fourier  transform 
arrangement  where  the  amplifier  crystal  of  Ba(NC)3)2  is  at  the  Fourier  transform  plane  of 
the  input  image. 
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A  standard  Air  Force  tri-bar  resolution  target  was  transilluminated  with  the  expanded 
1.56  pm  laser  pulses  from  the  Seed  Stokes  source  laser.  The  iUummatmg  laser  pulses 
were  greatly  attenuated  with  neutral  density  filters  up  to  a  total  optical  density  value  of 
9.8.  The  amplifier  pump  laser  was  timed  to  overlap  the  illuminating  pulses  inside  the 
Raman  amplifier  crystal.  The  output  energy  per  pulse  of  the  amplifier  laser  could  be 
adjusted  to  moderate  the  amount  of  gain  applied  to  the  illuminating  pulses.  Figure 
5.3. 1-2,  on  the  left,  shows  the  image  produced  when  no  amplifier  pump  pulses  were 
applied  to  the  SSRIA  crystal  (the  actual  output  gray  level  of  this  image  has  been  boosted 
so  that  the  tri-bars  can  be  seen  -  the  image  directly  out  of  the  camera  is  approximately  1- 
bit  over  background).  The  right  side  of  Figure  5. 3. 1-2  shows  the  corresponding  image 
with  the  amplifier  pump  laser  applied  to  the  SSRIA  crystal.  The  gray  level  is  now  at  the 
saturation  level  of  the  IR  camera  indicating  an  optical  gain  of  greater  than  500.  Also  note 
that  the  image  resolution  is  not  degraded  by  the  SSRIA. 
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Raw  Image  Amplified  Image 

(Full-Scale) 

Figure  5.3.1-2  Cpmparison  of  Unamplified  and  Raman  Amplified  Images  of  Air 

Force  Resolution  Target 

5.3.2  Range-gating 

The  tri-bar  resolution  target,  dluminated  by  the  seed  Stokes  beam,  is  imaged  into  the 
nonlinear  crystal.  The  seed  Stokes  image  light  is  amplified  by  stimulated  Raman 
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scattering  in  the  crystal.  Since  the  seed  Stokes  light  is  only  amplified  if  the  pump  light  is 
present  inside  the  crystal,  by  controlling  the  relative  timing  between  the  outgoing  seed 
Stokes  pulses  and  the  Q-switched  pump  pulses,  the  range  of  the  scene  imaged  can  be 
selected.  The  effective  “range-gate”  width  is  controlled  by  the  temporal  pulse  width  of 
the  pumps,  shows  the  range-gating  feature  of  SSRIA.  Starting  at  the  top  left  of  the 
figure,  the  upper  block  shows  a  single  tri-bar  image  that  is  unamplified  because,  as  is 
shovm  in  the  block  directly  below  it,  the  amplifier  pump  pulse  (the  wider  trace)  does  not 
overlap  the  illuminating  Stokes  pulse  (the  narrower  trace).  The  successive  images  and 
oscilloscope  traces  follow  to  the  right  and  then  fi:om  left  to  right  at  the  bottom.  These 
images  show  the  amplification  as  a  function  of  the  temporal  overlap  between  the 
illuminating  Stokes  pulses  and  the  amplifier  pump  pulses  in  the  Raman  crystal.  Each  set 
shows  the  result  of  shifting  the  timing  between  the  Stokes  pulse  and  the  amplifier  pump 
pulse  by  5  nanoseconds.  Note  the  mismatch  between  the  widths  of  the  Stokes  pulses  and 
the  amplifier  pump  pulses.  The  Stokes  pulses  are  nominally  5  nanoseconds  FWHM, 
while  the  amplifier  pump  pulses  are  nominally  30  nanoseconds  FWHM.  The  effective 
temporal  width  of  amplification  is  seen  to  be  about  15  nsec,  which  corresponds  to  a  range 
gatewidth  of  2.25  meters  in  air. 
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Figure  5.3.2-1  Range-gated  images  acquired  using  SSRIA. 

The  top  figures  show  the  amplified  images  of  an  AF  resolution  target  The  bottom 
figures  show  the  relative  transmitter  (red)  and  receiver  (green)  pulse  profiles  in  the 
plane  of  the  SSRIA.  The  relative  time  delays  are  denoted  in  5  ns  increments. 

5.3.3  Resolution 

For  the  experimental  demonstration  of  SSRIA  described  in  this  dissertation,  the 
system  MTF  is  a  product  of  the  MTF’s  due  to  diffraction,  optical  aberrations,  and  the 
SSRIA.  The  demonstration  was  designed  such  that  the  SSRIA  MTF  would  be  the 
limiting  term  in  the  optical  resolution.  Figtire  5.3.3-1  shows  a  theoretical  plot  of  the 
contribution  of  each  MTF  term  and  the  combined  MTF  for  the  breadboard  demonstration. 

In  the  experiment,  one  inch  optics  were  chosen  for  convenience.  The  clear  aperture  of 
the  Fourier  transform  lens  was  19.2  mm.  However,  since  the  Ba(N03)2  crystal  used  had 
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dimensions  of  10mm  x  10mm  x  44  mm,  the  limiting  aperture  stop  was  the  crystal  itself. 
The  optical  cutoff  frequency  for  a  limiting  aperture  of  1 0  mm  is : 


=  42.7  cycles  per  millimeter 


(5.2) 


The  MTF  contribution  due  to  optical  aberrations  was  estimated  by  a  raytrace  analysis 
of  the  SSRIA  optical  system.  The  geometrical  blur  size  at  the  IR  camera  was  determined 
to  be  approximately  8.6  |jm  in  diameter.  Assuming  that  the  optical  aberrations  result  in  a 
Gaussian  MTF  with  a  standard  deviation  of  a  =  05lr ,  where  r  is  the  blur  radius,  the 
MTF  due  to  optical  aberrations  is 

{5.3) 

whCTe  p  is  the  optical  frequ®cy  in  cycles  per  millimets’. 


Figure  5.3.3-1  Combined  MTF  of  Experimental  Demonstration 
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The  cutoff  frequency  of  the  SSRIA  is  controlled  by  the  diameter  of  the  pumped 
region,  d,  and  the  focal  length  of  the  inverse  Fourier  transform  \&as,f2,  by  the  equation: 


P  SSRIA  ~ 


f2^ 


(5.4) 


In  this  demonstration,  the  focal  length  was  chosen  for  convenience  and  the  pumped 
region  diameter  was  limited  by  the  pulse  energy  available  from  the  pump  laser.  For 
follow-on  work,  the  pumped  region  will  be  increased  to  5  mm  and  the  inverse  transform 
lens  win  be  chosen  to  demonstrate  a  cutoff  frequency  of  greater  than  40  line  pairs  per 
millimeter.  This  wiU  exceed  the  best  resolution  of  current  image  intensifier  tubes. 


Figure  5.3.3-2  Amplified  Image  of  Resolution  Group  +3 


Figure  5.3.3-2  shows  a  typical  amplified  image  of  the  resolution  target  with  Group  +3 
centered  in  the  picture.  In  this  image.  Element  5  was  resolvable.  Since  Group  +3, 
Element  5  corresponds  to  a  resolution  value  of  12.7  cycles  per  millimeter,  this  confirms 
the  MTF  model  for  the  SSRIA. 
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6  Raman  Image  Amplifier  System  Applications 

6.1  Eye-safe  Imaging  Lidar 

The  SSRIA  has  the  potential  to  improve  on  the  current  state  of  image  intensifier  and 
sti'eak  tube  perfonnance,  particularly  in  the  iofi’ared  where  photocathode  intensifiers  have 
poor  quantum  efficiency.  It  offers  the  potential  of  being  less  expensive,  simpler,  more 
rugged,  and  having  higher  performance  than  current  sensor  devices . 

Applications  for  military  use  include  eye-safe  imaging  lidars  that  can  be  used  for 
autonomous  vehicle  identification  and  targeting  at  20  km  and  greater  ranges,  such  as 
required  by  the  Air  Force’s  ERASER  program.  Reduced  background  radiation  through 
narrow  receiver  spectral  bandwidths  and  increased  range  resolution  from  short  amplifier 
laser  pulsewidths  are  realizable  with  SSRIA  seiisors.  The  increased  range  resolution  will 
also  improve  clutter  rejection.  SSRIA  also  has  the  potential  to  perform  optical  image 
processing  fimctions  that  will  allow  automatic  target  recognition  or  aimpoint  designation 
in  real-time  without  the  need  for  extensive  digital  signal  computations.  This  is  readily 
accomplished  in  SSRIA  by  pre-filtering  the  optical  image  before  amplification. 
Matched-filtered  output  fi:om  the  image  amplifier  can  be  obtained.  A  suitable  threshold 
can  then  be  applied  to  this  image  to  produce  an  all-optical  automatic  target  recognition 
process. 

6.2  Optical  Image  Processing 

In  most  lidar  system  today,  image  enhancement  and  automatic  target  detection  (ATD) 
are  performed  by  an  image  processing  computer  after  the  image  data  has  been  collected. 
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SSRIA  has  the  flexibility  to  perform  real-time,  Fourier-plane  pre-processing  of  optical 
images,  Figure  6.2-1  shows  a  schematic  diagram  of  how  this  can  be  accomplished. 
Active  and/or  passive  filters  can  be  added  in  a  Fourier-plane  of  the  image  to  filter  out 
undesirable  spatial  frequencies.  The  resulting  filtered  image  is  then  amplified  through 
the  SSRIA  and  collected  on  a  detector.  The  filters  can  be  passive  or  active  matched- 
filters  providing  a  real-time,  optical  processing  ATD  routine.  Phase-screens  could  also 
be  used  to  amplify  phase-objects  which  may  be  invisible  to  direct  imaging. 

Figure  6.2-1  shows  an  incoming  low-level  image  of  a  resolution  target.  The  optical 
filter  function  is  a  cross-shaped  bandpass  filter  with  a  central  ring  obscuration  to  filter  out 
the  low  spatial  frequency  intensity  variation  across  the  image.  This  filter  is  applied  in  an 
intermediate  Fourier-plaiie  of  the  image.  The  resulting  filtered  image  is  then  amplified 
by  SSRIA  and  the  final  amplified  image  shows  much  improved  contrast  of  the  resolution 


bars. 
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Figure  6.2-1  Filtered  Fourier-Plane  Image  Amplification 


Another  application  of  optical  image  processing  has  been  named  RAPOD  for  “Raman 
Amplified  Phase  Object  Detection.”  Certain  objects  are  invisible  to  standard  i-magiTig 
lidar  systems  because  they  involve  only  a  shift  in  optical  phase  rather  than  an  amplitude 
contrast.  These  objects  are  called  phase  objects  and  they  are  typically  associated  with  a 
varying  index  of  refi’action  over  the  imaged  scene.  In  air,  these  refractive  index 
variations  can  be  caused  by  a  sharp  pressure  shock  wave  such  as  the  supersonic  shock 
wave  caused  by  the  passing  of  a  supersonic  airplane  or  missile.  Other  phase  objects  of 
interest  are  strong  air  turbulence  and  the  wmgtip  vortices  shed  by  all  aircraft  wings  when 
generating  lift. 

Transparent  phase  objects  can  be  characterized  by  an  amplitude  transmittance  function 

t{x,y)  =  exp[i(p{x,y)]  (6.1) 
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where  (p(x,y)  is  the  spatially  varying  phase  shift  of  the  phase  object.  If  the  phase  shift  is 
smaller  than  1  radian,  the  amplitude  transmittance  can  be  approximated  by 

t{x,y)^\+i<p{x,y)  (6.2) 

where  terms  of  the  order  of  ^  and  higher  are  neglected.  The  first  term  of  Equation  (6.2) 
leads  to  a  strong  wave  component  that  passes  through  the  phase  object  witbout  change. 
The  second  term  generates  a  weaker  component  of  diffracted  light  that  is  deflected  away 
from  the  axis  of  the  system.  In  a  conventional  imaging  system,  the  image  produced  by 
the  phase  object  is  given  by 

/.=|l+/>f-l  (6.3) 

where  the  term  has  been  dropped  to  mamtain  consistency. 

In  1935,  Fritz  Zermke^'*  realized  that  the  diffracted  component  is  not  visible  because  it 
is  in  phase  quadrature  with  the  strong  background.  He  proposed  modifying  this  phase 
quadrature  relationship  by  mserting  a  phase-changing  plate  in  the  focal  plane.  The  strong 
background  component  is  brought  to  a  focus  on-axis  in  the  focal  plane  while  the 
diffracted  light  is  spread  away  from  the  focal  point.  By  placing  a  phase-changing  plate  in 
the  focal  plane  that  consists  of  a  glass  substr  ate  with  a  small  dielectric  dot  deposited  at 
the  focal  point  (i.e.  on  the  optical  axis),  the  phase  of  the  strong  background  component 
can  be  shifted  relative  to  the  diffracted  component.  For  example,  if  the  dot  causes  a 
phase  shift  of  7i/2  radians  relative  to  the  phase  retardation  of  the  diffracted  hght,  the 
intensity  in  the  image  plane  becomes 
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2 

=  ]f{\  +  <pf  ^l+2(p  (6.4) 

Therefore,  the  image  intensity  is  now  linearly  related  to  the  phase  shift  (p  and  exhibits 
positive  phase  contrast.  This  is  the  principle  of  the  Zemike  phase-contrast  microscope. 

Figure  6.2-2  shows  an  expanded  view  of  the  BAPOD  system.  The  RAPOD  system 
consists  of  a  crystal  of  Ba(N03)2  that  is  pumped  by  a  Nd:YAP  or  Nd:YAG  power 
oscillator  configuration.  The  phase-object  scene  illuminated  by  the  seed  Stokes  beam  is 
imaged  to  a  phase  contrast  plate  (located  ia  a  Fomier  plane),  which  has  a  small  dielectric 
dot  on-axis  for  shifting  the  phase  of  the  backgrormd  light  relative  to  the  diffracted  light  of 
the  phase  object.  The  phase-shifted  image  is  thea  relayed  to  a  second  Fourier  plane 
located  at  the  center  of  the  Raman  crystal.  The  seed  Stokes  image  tight  is  amplified  by 
stimulated  Raman  scattering  in  the  crystal.  Since  the  seed  Stokes  tight  is  oiily  amplified 
if  the  pump  tight  is  present  inside  the  crystal,  by  controlling  the  relative  timing  between 
the  outgoing  seed  Stokes  pulses  and  the  Q-switched  pump  pulses,  the  range  of  the  scene 
imaged  can  be  selected.  The  effective  “range  gate”  width  is  controlled  by  the  temporal 
pulse  width  of  the  pumps.  The  proposed  Q-switched  system  can  produce  pump  pulse 
widths  on  the  order  of  10  nanoseconds,  giving  an  effective  range-gate  on  the  order  of  1.5 
meters.  By  using  mode  locking  of  the  seed  Stokes  source  and  the  Nd:YAP  pump  laser,  it 
is  possible  to  produce  pulses  on  the  order  of  20-30  picoseconds,  resulting  in  effective 
range  gates  of  3-4  millimeters.  With  range  gates  this  short,  it  wiU  be  possible  to  image 
through  very  turbid  media,  amplifying  orily  the  unscattered  “ballistic”  photons  and 
rejecting  the  photons  which  have  been  multiply  scattered  by  the  media. 
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Figure  6.2-2  Raman  amplified  phase  object  detection  layout 


6.3  Marine  Lidar 

Other  applications  of  Raman  image  amplification  include  underwater  imaging  of 
mines  and  submarines.  Mine  countermeasure  (MCM)  systems  require  high  resolution, 
high  sensitivity  sensor  subsystems  to  perform  the  functions  of  mine  classification  and 
mine  identification.  LE)AR-based  miue  detection  systems,  which  are  currently  bemg 
added  to  the  Navy’s  inventory,  are  limited  in  their  ability  to  provide  identification  of  the 
mine-hke  objects  which  they  detect.  The  term  “identification”  in  this  context  means 
associating  a  sensed  mme-like  object  with  a  label  that  is  unique  to  the  target’s  external 
shape  and  configuration.  In  practical  terms,  this  means  that  the  sensing  element  must 
provide  sufficiently  high  spatial  resolution  so  that  the  image  of  a  particular  bottom  mine, 
for  example,  can  be  distinguished  from  a  discarded  oil  drum  or  some  other  bottom  debris. 

Range-gated  marine  LIDAR  systems  typically  employ  a  pulsed  laser  transmitter 
(usually  a  frequency  doubled  Nd:YAG  laser  at  a  wavelength  of  532  nm)  and  either 
intensified  CCD  cameras  or  photomultiplier  tubes  as  receivers.  Range  gating  is 
accomplished  by  controlling  the  timing  between  the  outgoing  laser  pulses  and  the  gating 
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or  exposure  of  the  sensing  element  of  the  LIDAR  system.  Image  intensifier  tubes  operate 
by  using  a  photoanissive  detector  coupled  to  a  microchannel  plate  (MCP)  electron 
amplifier.  Current  commercial  image  intensifiers  are  limited  to  gains  of  aroimd  10^, 
gating  times  of  aroimd  10  to  20  nanoseconds,  and  resolutions  of  around  50  to  100  pni  at 
the  photocathode.  The  resolution  is  primarily  limited  by  the  MCP  blur.  The  signal-to- 
noise  ratio  (SNR)  of  the  range-gated  image  is  degraded  by  the  quantum  efficiency  of  the 
photocathode  (typically  <  40%)  and  the  noise  figure  of  the  MCP  (a  degradation  of  SNR 
due  to  the  electron  multiplication  statistics  of  the  MCP). 

The  technique  of  SSRIA  can  outperform  current  image  intensifiers  and  offers  potential 
gains  of  10®,  gating  times  of  1-5  nanoseconds  (or  even  as  short  as  10-30  picoseconds  by 
using  laser  mode-locking  techniques),  and  resolution  of  5-10  pm  or  smaller.  The 
application  of  SSRIA.  to  the  marine  LIDAR  problem  is  called  MArine  Raman  Image 
Amplification  (MARIA). 

Figure  6.3-1  shows  a  schematic  of  the  MARIA  LIDAR  concept.  The  blue/green  laser 
transmitter  is  an  all-solid-state  laser  consisting  of  a  pumped  laser  crystal  and  an  external 
cavity,  second  harmonic  generator  (SHG).  The  optical  frequency  is  then  down-shifted  to 
the  first  Stokes  output  firequency  by  an  external  cavity  Raman  crystal.  The  Stokes  output 
illmninates  the  object  and  is  diffusely  reflected  back  to  the  receiver  telescope  optics.  The 
resulting  image  is  relayed  to  the  solid-state  Raman  image  amplifier  which  is  configured 
in  a  Fourier  transform  arrangement  where  the  amplifier  crystal  is  at  the  Fourier  transform 
plane  of  the  input  image. 


110 


Stokes  I  1. . t . I 


Input  *  •  ’ 

Image  Plane  Fourier  Transform  Image  Plane 

Plane 

Figure  6.3-1  MARIA  LIDAR  Concept 

The  Raman  image  amplifier  crystal  is  pvunped  by  a  beam  from  a  second  Q-switched, 
doubled  solid-state  laser.  The  Raman  amplification  process  only  occurs  when  the  pump 
light  is  present  in  the  Raman  crystal.  Therefore,  the  temporal  pulsewidth  of  the  amplifier 
pump  laser  determines  the  range-gatewidth  of  the  resulting  amplified  image. 

In  order  to  quantify  the  performance  difference  of  using  a  MARIA  LEDAR  system 
over  other  currently  available  underwater  imaging  systems,  a  set  of  common  parameters 
was  used  to  calculate  the  matched-filter  SNR  of  a  particular  target.  Both  1  nsec  and  3 
nsec  effective  gatewidths  were  evaluated  for  the  MARIA  system.  These  were  compared 
with  an  imaging  system  based  on  a  streak  tube  (STIL),  a  system  based  on  a  range-gated 
ICCD,  and  a  laser  line  scan  system.  The  LIDAR  system  parameters  used  for  MARIA  are 


shown  in  Table  6-1 . 


Table  6-1  MARIA  system  parameters 


Parameter 

Value 

Optical  wavelength,  X 

563  nm 

Laser  pulse  energy,  Q 

O.IJ 

Transmitter  half  angles,  0x,y 

0.011  rad 

Receiver  half  angles  (pixel),  0x,y 

0.0005  rad 

Receiver  area,  A^vr 

l^;r(0.18)" 

Amplifier  Noise  Factor  (see  Section 
2.4.2) 

2 

Raman  gain 

1000 

The  target  parameters  used  for  all  cases  is  given  in  Table  6-2. 


Table  6-2  Tar 

get  Parameters 

Parameter 

Value 

Radius,  po 

0.5  m 

Area,  Argt  =  npo^ 

0.2571 

Reflectance,  R 

0.08 

The  parameters  used  for  the  STIL  are  given  in  Table  6-3. 


Table  6-3  STIL  system  parameters 


Parameter 

Value 

Optical  wavelength,  X 

532  nm 

Laser  pulse  energy,  Q 

0.1  J 

Transmitter  half  angles. 

Along  track:  0^ 

0.0005  rad 

Cross  track:  0y 

0.25  rad 

Receiver  half  angles  (pixel), 

Along  track:  0* 

0.0005  rad 

Cross  track:  0y 

0.0005  rad 

Receiver  area,  A^cvr 

y^7c{0X%f 

Photodetector  quantum  efficiency,  T|q 

0.1 

The  systeni  parameters  for  the  ICCD  system  are  given  in  Table  6-4. 


Table  6-4  ICCD  system  parameters 


Parameter 

Value 

Optical  wavelength,  X 

532  nm 

Laser  pulse  energy,  Q 

0.1  J 
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Transmitto’  half  angles,  0x,y 

0.011  rad 

Receiver  half  angles  (pixel),  6x,y 

0.0005  rad 

Receiver  area,  Arcvr 

^;r(0.18)^ 

Photodetector  quantum  efficiency,  T]q 

0.1 

Finally,  the  system  parameters  used  for  the  laser  line  scan  case  are  given  Table  6-5. 


Table  6-5  Laser  line  scan  system  parameters 


Parameter 

Value 

Optical  wavelength,  X 

532  nm 

Laser  power,  P 

1.5  W 

Transmitter  half  angles,  9x,y 

0.00025  rad 

Transmitter  offset,  yxmtr 

-1.0  m 

Receiver  half  angles  (pixel),  0x,y 

0.032  rad 

Receiver  offset,  Yrevr 

1.0  m 

Receiver  area,  A^vr 

m^ 

Photodetector  quantum  efficiency, 

0.1 

Full  scan  angle,  @ 

0.5  rad  (28  deg) 

Scan  time 

0.0667  sec 

AH  system  pai'ameters  were  chosen  to  be  representative  of  current  state-of-the-art 
devices.  The  scan  time  and  laser  power  of  the  laser  line  scan  system  was  selected  to 
deposit  the  same  laser  energy  over  the  scan  as  used  with  the  other  systems.  This  means 
that  any  differences  in  SNR  will  be  due  to  the  detailed  detection  processes  of  each  system 
and  not  due  to  any  differences  in  Hlummation. 

Figure  6.3-2  shows  the  Jerlov^^  ocean  water  types  with  1, 13,  and  m  being  clear  ocean 
water  types  and  1,  4,  5,  7,  and  9  being  coastal  water  types.  Notice  that  the  irradiance 
transmittance  for  Jerlov  coastal  type  4  is  fairly  flat  between  532  nm  and  563  nm.  Thus 
there  is  no  significant  difference  between  using  the  Stokes-shifted  wavelength  at  563  nm 
and  the  standard  doubled  Nd;YAG  wavelength  at  532  nm  for  Jerlov  coastal  type  4.  This 
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is  trae  for  most  of  the  coastal  water  types,  with  the  more  turbid  cases  such  as  types  7  and 
9  being  slightly  more  transmissive  at  563  nm. 


0 


Doubled  Nd:YAG 

Doubled  Nd:YAG 
and  Ba(N03)2 


Figure  6.3-3  shows  the  performance  comparison  between  all  the  system  types.  The 
differences  between  the  various  curves  are  primarily  due  to  two  factors.  The  first  factor 
is  the  effective  gatewidth  employed  by  the  system.  The  jBrst  curve,  MARIA  (1  nsec), 
represents  the  optimum  performance  which  can  be  obtained  in  this  type  of  coastal  water. 
The  temporal  spreading  due  to  multiple  scattering  is  just  under  1  nsec  for  the  round  trip 
through  the  water  column.  Thus  any  shorter  effective  gatewidth  would  r^ult  in  a  lower 
SNR.  To  achieve  this  performance,  the  MARIA  system  would  have  to  employ  pulses  of 
around  1  nsec  FWHM.  This  would  present  a  major  challenge  to  the  aniplifier  pulse  laser. 
The  second  curve,  MARIA  (3  nsec),  would  employ  pulses  of  less  than  3  nsec  FWHM. 
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This  is  still  challengiag  for  the  amplifier  pulse  laser  but  can  be  achieved  through  the  use 
of  cavity-dumping  techniques. 

The  difference  between  the  MARIA  (3  nsec)  and  the  STIL  (3  nsec)  performance  is 
due  to  the  second  factor,  which  is  the  effective  quantum  efficiency  of  detection.  As 
discussed  in  the  Amplifier  Theory  Section,  the  MARIA  system  has  an  overall  effective 
quantum  efficiency  of  50%.  That  is,  the  SNR  noise  power  factor  approaches  the  limiting 
value  of  2  as  the  gain  is  increased.  In  the  STTL  system  the  effective  quantum  efficiency 
is  10%  due  to  the  lower  quantum  efficiency  of  the  photocathode  and  the  noise  factor  of 
theMCP. 

The  effective  range  gate  used  for  the  ICCD  system  is  25  nsec  (or  about  9  feet).  This 
assumes  a  transmitter  pulsewidth  of  15  nsec,  which  is  typical  for  the  underwater  imaging 
LEDAR  systems  in  use  today 

Given  a  criterion  for  detection  of  SNR  =10  dB,  the  target  detection  depths  achievable 
with  each  system  are  given  in  T able  6-6 . 


.  115 


-H- MARIA  (1  nsec) 
MARIA  (3  nsec) 
-^STIL  (3  nsec) 
-•-ICCD  (25  nsec) 
Line  Scan 


Figure  6.3-3  Performance  Comparison  of  MARIA  with  STBL,  ICCD,  and  Laser 

Line  Scan  Systems 


Table  6-6  Target  Detection  Depths  (SNR  =  10  dB  criterion) 


System 

Detection  depth,  D 

Diffusion  lengths,  KD 

MARIA  (1  nsec) 

37  m 

8.0 

MARIA  (3  nsec) 

36 

7.8 

STIL  (3  nsec) 

32 

6.9 

ICCD 

30 

6.5 

Laser  line  scan 

23 

5.0 

A  design  study  was  performed  to  determine  the  optimum  design  for  the  Stokes  seed 
laser  used  to  illuminate  underwater  targets.  Four  candidate  pulsed  laser  designs  were 
considered:  1.  Flashlamp-pumped,  Q-switched  Nd.YAG;  2.  Flashlamp-pumped,  Q- 
switched  and  cavity  dumped  Nd:YAG;  3.  Flashlamp-pumped,  Q-switched  and  mode 
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locked  Nd:YAG;  and  4.  Flashlamp-piunped,  Nd:YAG  regenerative  amplifier  (regen). 
RasWamp-pumped  lasers  were  chosen  for  their  low  relative  cost. 

The  standard  Q-switched  Nd:YAG  laser  design  is  shown  in  Figure  6.3-4.  All  of  the 
Stokes  seed  laser  designs  are  shown  with  the  second  harmonic  generator  (SHG)  and  the 
external  Raman  crystal.  The  pulse  shape  of  the  Q-switched  Nd;YAG  is  governed  by  the 
gain  characteristics  of  the  laser  medium  and  the  pulse  width  is  typically  10  -  20  nsec. 
With  high  optical  pumping  and  by  shortening  the  resonator  as  much  as  possible,  the 
lower  limit  of  a  standard  Q-switched  Nd:  YAG  is  about  8  nsec.  This  temporal  limit  is  far 
greater  than  the  pulse  stretching  caused  by  multiple  scattering  in  typical  ocean  waters. 
Thus  the  full  SNR  potential  of  range-gating  is  not  reached  by  a  standard  Q-switched 
Nd:YAG. 


PC  p 


HR@1.06pm 


Nd:YAG 


SHG 


OC@1.06  pm 

.  Figure  6.3-4  Q-Switched  Nd:YAG 


Figure  6.3-5  shows  the  second  candidate  laser  design,  the  Q-switched  and  cavity- 
dumped  Nd;YAG.  In  a  cavity  dumped  laser,  both  end  mirrors  of  the  laser  cavity  are 
essentially  100%  reflectors.  When  the  flashlamp  is  fired,  the  Q-s witch  is  held  at  a  low 
retardation  voltage,  which  results  in  the  hoiizontally  polarized  fluorescence  being 
transmitted  through  the  thin-fihn  polarizer  and  prevents  lasing.  Upon  reaching  peak 
energy  storage  in  the  Nd:YAG  rod,  the  Pockels  cell  is  biased  to  its  half-wave  retardation 
voltage.  The  resulting  vertically  polarized  Hght  is  now  reflected  by  the  polarizer  to  the 
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other  HR  and  the  laser  pulse  energy  begins  to  build  up.  When  the  Q-switched  pulse 
power  in  the  cavity  reaches  a  maximum,  the  bias  is  quickly  removed  from  the  Pockels 
cell.  The  cavity  energy  is  then  transmitted  through  the  polarizer,  essentially  dumping  the 
entire  cavity  within  one  round-trip.  The  width  of  the  laser  pulse,  tp,  is  then  given  by: 


(6.5) 


c 

and  for  a  short  0.5  m  cavity  is  about  3.3  nsec.  This  cavity-dumped  design  can  produce 
almost  as  much  energy  per  pulse  as  the  standard  Q-switched  design  and  is  operated  at  the 
same  range  of  PRF.  The  cost  is  somewhat  higher  due  to  the  more  demandmg  Q-switch 
electronics  required  to  switch  quickly. 


PC  X/4 

Polarizer 

HR@1.06|im 

Figure  6.3-5  Q-Switched  &  Cavity  Dumped  NdrYAG 
Figure  6.3-6  shows  the  third  candidate  design,  the  Q-switched  and  mode  locked 
NdiYAG.  By  placiag  an  amplitude  modulator  inside  the  laser  cavity  and  diivmg  it  at 
exactly  the  frequency  separation  of  the  axial  modes,  the  laser  can  be  induced  to  generate 
a  train  of  mode-locked  pulses  with  a  PRF  of  f^^cjlL.  The  inclusion  of  a  Q-switch  in 
the  laser  cavity  results  m  a  burst  of  mode-locked  pulses  that  are  contained  within  the 
envelope  of  a  100  -  200  nsec  Q-switch  pulse.  The  individual  mode-locked  pulses  are 
between  10  and  30  psec  FWHM  with  an  energy  of  up  to  10  mJ  at  the  peak  of  the  Q- 
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switch  envelope.  These  nairow  pulses  will  be  rapidly  spread  by  the  multiple  scattering  in 
ocean  water  and  will  approach  1  nsec  in  width  on  the  return  path.  Thus  the  temporal 
(and,  therefore,  range)  resolution  afforded  by  using  30  psec  pulses  cannot  be  realized  in 
an  underwater  environment.  In  addition,  the  PRF  of  the  mode-locked  pulses  is  much  too 
frequent  to  allow  pulses  to  propagate  to  an  underwater  target  and  return  before  the 
emission  of  another  pulse.  Thus  the  returns  of  many  pulses  at  many  different  ranges  wiU 
be  superimposed,  causing  SNR  degradation. 


PC  p 


Nd:YAG 


AO 

Modulator 


SHG  Raman 


HR@1.06|mi  OC@1.06  imi 

Figure  6.3-6  Q-Switched  &  Mode  Locked  Nd:YAG 
The  fourth  candidate  design,  the  regenerative  amplifier,  allows  a  way  to  produce 
individual  picosecond  pulses  at  lower  PRF’s,  thus  avoiding  the  overlap  of  pulse  returns  at 
different  ranges.  This  laser  system  consists  of  a  mode-locked  oscillator,  which  provides 
low-energy,  mode-locked  pulses,  and  a  regenerative  amplifier,  which  selects  and  traps  an 
individual  mode-locked  pulse,  amplifies  it  greatly,  and  switches  it  back  out.  To 
understand  the  operation  of  the  regenerative  amplifier,  consider  a  horizontally  polarized 
pulse  from  a  mode-locked  oscillator  that  is  transmitted  through  the  polarizer  in  Figure 
6.3-7.  The  horizontally  polarized  pulse  wUl  experience  a  45°  rotation  at  the  X/2  plate  and 
a  further  45°  rotation  from  the  Faraday  rotator.  The  now  vertically  polarized  pulse 
reflects  from  the  polarizer  into  the  laser  resonator.  After  passing  through  the  X/4  plate 
and  reflecting  from  the  HR,  the  pulse  once  again  passes  through  the  X/4  plate  resulting  in 
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a  horizontally  polarized  pTilse  traveling  toward  the  right.  As  soon  as  the  pulse  leaves  the 
Pockels  cell  traveling  toward  the  right,  the  Pockels  cell  voltage  is  quickly  switched  to  X/4 
retardation.  Therefore,  when  the  pulse  next  makes  its  round  trip  through  the  Xy4  plate 
and  Pockels  cell,  it  remains  horizontally  polarized  and  is  trapped  inside  the  cavity.  After 
many  passes  in  the  resonator,  the  energy  of  the  pulse  reaches  a  maximum  due  to  gain 
saturation.  At  this  point,  the  Pockels  cell  is  switched  to  yjl  retardation  voltage.  The 
pulse  now  experiences  a  90°  polarization  rotation  and  the  vertically  polarized  pulse  is 
reflected  by  the  two  polarizers  to  the  output  port  of  the  laser  system. 

The  regenerative  amplifier  design  can  produce  30  psec  pulses  at  a  repetition  rate  up  to 
10  kHz  with  a  pulse  energy  of  up  to  5  mJ.  Although  the  regen  solves  the  problem  of 
overlapping  mode-locked  pulses  as  exhibited  by  the  third  design,  the  pulse  energy  is  still 
small  and  the  30  psec  pulses  are  rapidly  stretched  to  1  nsec  by  the  ocean  water. 


PC  V4  Polarizer  Nd:YAG  HR@1.06nm 


Mode  Locked  Oscillator 

Figure  6.3-7  Regenerative  Amplifier 

Table  6-7  shows  a  summary  comparison  of  the  four  candidate  Stokes  seed  laser 
designs.  Based  on  the  desire  to  match  the  temporal  spreading  expected  from  multiple 
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scattering  in  the  ocean  and,  thereby,  to  realize  the  optimum  performance  in  terms  of 
range  resolution  and  water  backscatter  rejection,  the  best  candidate  for  MARIA  is  the  Q- 
switched  and  cavity  dumped  design.  This  design  will  provide  the  best  overall 
performance  within  reasonable  bounds  for  cost. 


Table  6-7  Laser  System  Trade-off  Comparison 


Laser  Pulsing 
Type 

Pulsewidth 

(nsec) 

PRF 

(Hz) 

Pulse  Energy 
(mJ) 

Relative 

Cost 

Q-switched 

8-10 

1-100 

100-500 

Low 

Q-switched  & 
cavity  dumped 

3-5 

1-100 

100-200 

Medium 

Q-switched  & 
mode  locked 

10  -  30  psec 

100-300 

MHz 

1-10 

Medium 

Regen 

10-30psec 

1-lOkHz 

<5 

Hgh 

The  first  Stokes  output  will  be  frequency  doubled  using  an  uncoated  8x8x10  mm^ 
LBO  crystal  cut  for  type  I  non-critical  phase  matching  (0=90°,  (1)=0°).  The  phase 
matching  temperatui'e  of  42°C  is  maintained  by  enclosing  the  crystal  in  a  temperature 
controlled  oven.  LBO  was  chosen  as  the  doubling  material  because  it  exhibits  a  superior 
acceptance  bandwidth  (AT-L  s  6°  C-cm)  and  damage  threshold  (>1  GW/cm2)  over  other 
candidate  materials.  Previous  experimental  results  have  shown  that  the  doubling 
efficiency  is  around  55%. 

The  cutoff  fi-equency  of  the  MARIA  is  controlled  by  the  diameter  of  the  pumped 
region,  d,  and  the  focal  length  of  the  inverse  Fourier  transform  lens,/2,  by  the  equation: 


P  MARIA  ~ 


d 


(6.6) 
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The  pumped  region  will  be  increased  to  5  mm  in  diameter  and  the  inverse  transform 
lens  will  be  chosen  to  demonstrate  a  cutoff  frequency  of  greater  than  40  line  pairs  per 
millimeter.  This  will  exceed  the  best  resolution  of  current  image  intensifier  tubes. 

The  limiting  resolution  of  the  receiver  is  only  half  of  the  story  for  a  marine  imaging 
system.  The  resolution  of  images  presented  to  the  receiver  is  degraded  by  other  factors 
related  to  the  propagation  of  light  in  the  turbid  medium  of  ocean  water.  First,  there  is  the 
image  blur  due  to  the  forward  scattering  nature  of  tbe  intervening  water.  This  effect  can 
be  quantified  by  the  point  spread  function  of  the  water  column  (also  called  the  beam 
spread  function).  Essentially,  tiie  reflected  light  distribution  from  the  illuminated 
object  is  convolved  with  the  point  spread  function  of  the  water,  thereby  blurring  the 
image.  This  phenomena  affects  any  type  of  imaging  system,  be  it  a  range-gated  imager 
or  a  laser  line  scan  system.  The  second  effect  of  imaging  through  water  is  that  water 
backscatters  the  illuminating  laser  beam  and  this  backscattered  light  adds  to  the  light 
reflected  from  a  target,  thereby  washing  out  the  inherent  contrast.  It  is  this  effect  which 
all  active  marine  imaging  systems  seek  to  minimize.  In  a  range-gated  system,  such  as 
MARIA,  the  backscattered  light  is  minimized  by  only  opening  the  receiver  to  receive  the 
light  reflected  from  a  definite  range.  By  minimizing  the  width  of  the  range  gate,  the 
amount  of  backscattered  light  is  minimized.  This  has  the  effect  of  minimizing  the  veiling 
light  due  to  the  surrounding  water,  and  maximizing  target  image  contrast.  In  a  laser  line 
scanned  system,  the  backscattered  light  is  minimized  by  using  a  bistatic  configuration. 
The  laser  output  aperture  is  physically  separated  from  the  receiver  aperture  and  their 
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respective  fields-of-view  only  overlap  near  the  target.  Any  water  appearing  in  the 
overlap  backscatters  light  to  the  receiver  and  degrades  contrast. 

A  recent  study  by  Ronald  E.  Walker  and  John  W.  McLean^^  has  dmved  lidar 
equations  for  turbid  media  with  the  inclusion  of  pulse  stretching  caused  by  multiple 
scattering.  Included  in  this  work  was  an  analytical  comparison  of  the  contrast  transfer 
function  (CTF)  of  a  range-gated  imaging  (RGI)  system  with  the  CTF  produced  by  a 
comparable  laser  line  scanned  (LLS)  system.  Figure  6.3-8  shows  the  results  of  this 
computation  for  the  Coastal  Ocean  case  presented  in 
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Table  6-8.  The  contrast  and  SNR  are  plotted  as  a  function  of  spatial  frequency  for  a 
depth,  D,  corresponding  to  4,  6,  and  8  beam  attenuation  lengths  for  the  CO  case  (i.e., 
c  •  £>  =  4, 6,  or  8  ).  As  shown  in  Figure  6 .3-8 ,  the  CTF  produced  by  both  the  imaging  Hdar 
system  and  the  laser  line  scan  system  has  a  cutoff  near  1  cycle/cm. 


Figure  6.3-8  SNR  and  Contrast  (CTF)  calculations  for  shot-noise-limited  LLS 
(solid)  and  RGl  (dashed)  systems  at  4,  6  and  8  beam  attenuation  lengths  in  Coastal 
Ocean  (CO)  waters. 
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Table  6-8  Summary  of  Water  Types 


Water  Type 

c  (m") 

COo 

k  (m") 

Coastal  Ocean  (CO) 

0.398 

0.551 

0.22 

Tuitid  Harbor  (TH) 

2.19 

0.833 

0.67 

(c  is  the  beam  attenuation  coefficient,  coo  is  the  single  scattering  albedo, 
and  k  is  the  diffuse  attenuation  coefficient  of  ocean  water) 


Figure  6.3-8  shows  that  at  4  beam  attenuation  lengths  there  is  very  little  difference  in 
the  two  types  of  imaging  systems.  However,  at  6  and  8  beam  attenuation  lengths  the 
range-gated  imaging  system  outperforms  the  LLS  system  in  terms  of  image  contrast. 
This  is  a  consequence  of  the  increased  overlap  of  the  transmitter  and  receiver  fields-of- 
view  in  the  LLS  at  longer  ranges.  The  increased  overlap  leads  to  more  backscattered 
light  admitted  to  the  image.  The  RGI  system  rejects  the  backscattered  light  through  use 
of  a  narrow  range-gate,  which  is  essentially  independent  of  the  object  range  (pulse 
stretching  by  the  water  may  cause  some  increase  of  backscatter  with  range). 

Figure  6.3-10  shows  the  SNR  and  CTF  comparison  in  Turbid  Harbor  (TH)  waters. 


TM0eft.C7FKK.GRF 


Figure  6.3-10  SNR  and  Contrast  (CTF)  calculations  for  shot-noise-limited  LLS 
(solid)  and  RGI  (dashed)  systems  at  4,  6  and  8  beam  attenuation  lengths  in  Turbid 
Harbor  (TH)  waters. 
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As  shown  in  Figure  6.3-10,  although  the  ultimate  spatial  frequency  cutoff  for  the  two 
systems  is  the  same  (aroimd  10  cycles/cm),  the  RGI  system  shows  much  reduced  contrast 
at  intermediate  spatial  frequencies  due  to  the  increased  amount  of  glow/blur  caused  by 
forward  scattering  of  light  emanating  from  surrounding  areas  of  the  target  being  imaged 
into  other  areas.  Because  of  the  narrow  field-of-view  of  the  LLS  system,  there  is  less 
glow  present  and,  therefore,  the  image  contrast  is  higher.  In  effect,  the  LLS  system  is 
equivalent  to  a  high-pass  filter  on  the  image.  The  forward  glow/blur  manifests  itself  as  a 
low  spatial  frequency  noise  source.  Forward  glow  rejection  can  be  implemented  in  LRG 
systems  by  affecting  the  same  high-pass  filter  b^dre  detection  as  the  LLS  systems.  In 
MARIA  this  can  be  done  by  using  an  optical  high-pass  filter  (see  Section  6.2).  The 
MLARIA  optical  architecture  utilizes  Fourier  Transform  (FT)  imaging,  which  allows  real¬ 
time  optical  pre-processing  of  the  signal  (image)  via  FT-plane  intensity  and/or  phase 
masks.  Standard  ICCD-based  LRG  systems  cannot  use  this  blur/glow  rejection  technique 
because  the  filtering  can  only  be  performed  after  detection,  and  therefore  after  detector 
noise  has  been  added.  Therefore,  MARIA  offers  the  blur/glow  rejection  capabilities  of  the 
LLS,  but  with  the  single-shot  imaging  capabilities  of  the  LRG. 

The  CTF  of  a  marine  imaging  system  depends  not  only  on  the  turbidity  of  the  ocean 
water,  but  also  on  the  specific  fields-of-view,  effective  gatewidths,  pulsewidths,  and  other 
system  parameters.  It  is  not  clear  that  either  the  imaging  lidar  or  the  laser  line  scan 
system  is  generally  the  winner  in  terms  of  contrast  performance  or  image  resolution.  It  is 
clear  that  a  marine  imaging  system  based  on  MARIA  can  produce  high  resolution  images 
at  competitive  ranges  (or  dejpths)  in  most  of  the  required  ocean  waters. 


6.4  Medical  Lidar 
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Medical  applications  would  include  imaging  through  highly  scattering  media.  The 
SSRIA  sensor  would  enable  amplification  and  detection  of  ballistic  photons  transmitted 
through  tissue  with  very  high  range  resolution.  This  resolution  would  enable  production 
of  3  dimensional  images  of  small  millimeter  dimension  targets  such  as  tumors. 

A  modification  of  Raman  imaging  technology  has  potential  for  imaging  sub-surface 
skin  carcinomas.  By  using  mode-locking  technology,  the  Raman  illuminator  and 
amplifier  pump  laser  pulses  can  be  temporally  compressed  to  tens  of  picoseconds 
FWHM.  The  resulting  gate  for  Raman  imaging  would  allow  separation  of  the  “ballistic” 
photons,  which  pass  unscattered  through  soft  tissue,  firom  the  multiply  scattered  photons, 
which  would  normally  mask  the  desired  image.  Mahon,  et  al.,  have  demonstrated  that 
features  as  small  as  125  |mi  can  be  imaged  through  dense  scattering  media  with  light 
extinction  factors  of  up  to  e^^  using  Raman  image  amplification  in  Hydrogen  gas.^^  By 
using  solid-state  Raman  image  amplification,  a  medical  imaging  device  can  be  made 
smaller  and  more  compact,  with  even  higher  spatial  resolution.  The  technology  transfer 
could  result  in  much  earlier  detections  of  carcmomas,  particularly  those  associated  with 


Melanomas. 


7  Conclusions 


In  conclusion,  this  dissertation  has  shown  that  Solid-state  Raman  Image  Amplification 
is  a  feasible  technology  for  producing  range-gated  imagery  in  the  eye-safe  infrared. 
SSRIA  is  capable  of  amplifying  low-level  optical  images  ,  with  gains  of  10^  to  10^  with 
the  addition  of  only  quantum-limited  noise.  The  high  gains  available  from  SSRIA  can 
compensate  for  low  quantum  efficiency  detectors  in  the  SWIR  and  can  reduce  the  need  to 
cool  detectors  to  reduce  thermal  noise.  The  range-gate  of  SSRIA  was  shown  to  be 
controlled  by  the  pulsewidth  of  the  amplifier  pump  laser  and  can  be  made  as  short  as  30  - 
100  cm,  using  pump  pulses  of  2-6.7  nsec  FWHM.  The  use  of  SSRIA  in  an  imaging 
lidar  system  has  been  shown  to  result  m  higher  SNR  images  throughout  a  broad  range  of 
incident  light  levels,  in  contrast  to  the  increasing  noise  factor  occurring  with  reduced  gain 
in  ICCDs.  General  equations  for  the  gain,  noise,  and  optical  design  of  SSRIA  have  been 
presented.  A  theoretical  framework  for  the  optical  resolution  in  SSRIA  has  been 
presented  and  it  was  shown  that  SSRIA  can  be  designed  to  produce  higher  resolution  than 
is  att^able  by  ICCDs.  SSRIA  was  also  shown  superior  to  ICCDs  in  rejecting  unwanted 
sunlight  background,  further  increasing  the  SNR  of  images.  Lastly,  it  was  shown  how 
SSRIA  can  be  combined  with  optical  pre-filtering  to  perform  optical  image  processing 
fimctions  such  as  image  enhancement  and  automatic  target  detection/recognition. 
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